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Abstract. We present the first general relativistic hydrodynamic models of the launch and evolution of relativistic 
jets and winds, driven by thermal energy deposition, possibly due to neutrino-antineutrino annihilation, in the 
close vicinity of black hole-accretion torus systems. The latter are considered to be the remnants of compact 
object mergers. Our two-dimensional simulations establish the link between models of such mergers and future 
observations of short gamma-ray bursts by the SWIFT satellite. They show that ultrarelativistic outflow with 
maximum terminal Lorentz factors around 1000 develops for polar energy deposition rates above some 10** ergs~^ 
per steradian, provided the merger environment has a sufficiently low baryon density. By the interaction with the 
dense accretion torus the ultrarelativistic outflow with Lorentz factors F above 100 is coUimated into a sharp-edged 
cone that is embedded laterally by a wind with steeply declining Lorentz factor. The typical semi-opening angles 
of the F > 100 cone are 5° — 10°, corresponding to about 0.4 — 1.5% of the hemisphere and apparent isotropized 
energies (kinetic plus internal) up to ^ 10^^ erg although at most 10 — 30% of the deposited energy is transferred 
to the outflow with F > 100. The viability of post-merger black hole-torus systems as engines of short, hard 
gamma-ray bursts is therefore confirmed. The annihilation of neutrino-antineutrino pairs radiated from the hot 
accretion torus appears as a suitable energy source for powerful axial outflow even if only « 10**^ erg are deposited 
within a cone of 45° half-opening angle around the system axis. Although the torus lifetimes are expected to be 
only between some 0.01 seconds and several 0.1 seconds, our models can explain the durations of all observed 
short gamma-ray bursts, because different propagation velocities of the front and rear ends will lead to a radial 
stretching of the ultrarelativistic flreball before transparency is reached. The ultrarelativistic flow reveals a highly 
non-uniform structure caused by the action of Kelvin-Helmholtz instabilities that originate at the fireball-torus 
interface. Large radial variations of the baryon density (up to several orders of magnitude) are uncorrelated with 
moderate variations of the Lorentz factor (factors of a few) and fiuctuations of the gently declining radiation- 
dominated pressure. In the angular direction the Lorentz factor reveals a nearly fiat plateau-like maximum with 
values of several hundreds, that can be located up to 7° off the symmetry axis, and a steep decrease to less than 
10 for polar angles larger than 15°-20°. Lateral expansion of the ultrarelativistic core of the fiow is prevented 
by a subsonic velocity component of about 0.05c towards the symmetry axis, whereas the moderately relativistic 
wings show a subsonic sideways infiation with less than 0.07c (measured in the frame comoving with the radial 
fiow). 
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1. Introduction 

Gamma-ray bursts (GRBs) are transient astrophysical 
events vi^hich are characterized by the release of enormous 
amounts of energy in non-thermal, highly variable (in time 
and event by event) emission on timescales between some 
milliseconds and about a thousand seconds. The distribu- 
tion of burst durations reveals two classes, "short bursts" 
which last less than about 2 seconds and have typical du- 
ration of 0.3-0.5 seconds, and "long bursts" with typically 
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softer spectra and a mean du rations of about 30-60 sec- 
onds teouveliotou et al.lll993ll . 

So far only long bursts could be sufficiently well local- 
ized on the sky to enable multi- wavelength observations. 
In many such cases afterglow emission was detected, the 
host galaxies were discovered, and the associated measure- 
ment of spectral line redshifts confirm ed the cosmologi- 
cal distances to t he burst sources (e. g. , 'Costa et al."l 997t 
Piro et alJ Il998t Ivan Paradiis et all 11997: Bremer etTll 
1998t iFrail eraLlll997l) . Observational evidence also indi- 
cates that the bursts originate from star formation regions, 
suggesting a possible link to massive star explosions. This 
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idea was nourished by several cases where a bump-hke 
re-brightening phase of the decaying afterglow emission 
was interpreted as a superimposed supernova lightcurve. 
The link was directly confirmed by two cases where the 
spatial and temporal coincidence of a GRB with a pecu- 
liar kind of sup ernova explosion wa s found, GRB980425 
with SN1998bw ("Ga lama etaDllQQSll . wh ich was still am- 
biguous, and G RB030329 with SN2003dh l)Stanek fc et all 
bnontlHiorth et alJl200^ . For short bursts such constrain- 
ing observations do not exist. 

According to the most widely accepted theoretical con- 
cept GRBs are understood as the consequence of an ex- 
plosive release of energy associated with the accretion of 
matter at huge rates (up to a several solar masses per sec- 
ond) onto stellar-mass black holes, e.g., as a consequence 
of the birth of such a black hole (BH) by the catastrophic 
collapse of a massive star in a so-called "collapsar" event 
llWooslevlll993t iMacFadyen fc Wooslelll999t IaIov et alJ 
l20n(i iMacFadvenet all I2OOII) or by a merg e r of com- 
pact objects in a close binarv llPacvhskilll98f : iGoodmanI 
11984 iFichler et a].1ll989t iMochkovitch et a,].lll99,-^ . The 
nascent black hole is initially girded by a thick gas torus 
from which it swallows matter at a hypercritical rate, 
i.e., at a rate far in excess of the Eddington limit. In 
this case radiation is advected inward with the accre- 
tion flow an d the cooling is dom inated by the emission 
of neutrinos ijPopham et al.lll999l) . a situation for which 
the term "n eutrino-dominated ac cretion disk" (NDAF) 
was coined l|Naravan et alJ l200lf) . The annihilation of 
neutrino-antineutrino pairs in the immediate vicinity of 
the black hole-torus system is considered to create an 
e"'"e~-pair plasma-photon fireball which could power an 
ultrarelativistic outflow of baryons (with typical Lorentz 
factors > 10 ^) provided the baryo n loading remains suffi- 
ciently low ( Gavallo fc R,eejll978fl . Alternatively, magne- 
tohydrodynamic processes might accelerate matter to ul- 
trarelativistic velocities by tapping the gravitatio nal bind- 
ing energy of the accretio n flow (e.g., IPrenkhahn .2002; 
iDrenkhahn fc Spruitll2002l) or the rot ational energy of the 
black hole l|Blandford fc Znaie'5ll977l) . When these highly 
relativistic ejecta finally reach optically thin conditions at 
much larger radii (~ 10^^ cm) the kinetic energy of the 
outflow is dissipated and the observable GRB is produced 
mainly by synchrotron radiation and i nverse Compto n 
scattering (e.g., see the recent review bv lMeszarojl2002|) . 

A large number of theoretical studies focused on the 
physics of the relativistic fireball and the formation of the 
observable properties of the GRB and of the afterglow 
emission (for reviews, see e.g.. lMeszarosl20o3|Piranll999|) . 
assuming that the energy release by some compact cen- 
tral engine has been able to produce the ultrarelativistic 
outflow. However, only few of the suggested astrophysical 
scenarios or concepts have so far been investigated by de- 
tailed numerical simulations which attempt to develop a 
quantitative picture of the physical processes that play a 
role at the source of the energy and that are responsible 
for driving the highly relativistic ejection of matter. 



Hypercritical steady-state accretion from accretion 
disks and tori was flrst discu ssed with semi-analytic 
mo dels bv IPopham et al. (l999 !), and later recons idered 
bv iNaravan et alJ ll2001.) .Pi Matteo et all l|2002^ . and 
iKohri fc Mineshigg l(2002) . These works have brought 
a theoretical understanding of the conditions that 
are present in the vicinity of the accreting BH and 
which determine the efficiency of energy loss by neu- 
trino emission and the efficiency of energy conver- 
sion by neutrino-antineutrino annihilation. To this end, 
three-dime nsional hydrodynamic simulations were per- 
formed by Ruffcrt fc J aiikal l)l999() and more recently by 
ISetiawan et al.. 1,2004;) for the time- dependent accretion in 
BH-tori systems as the remnants of neutron star-neutron 
star (NS-I-NS) and NS-I-BH mergers. In this case the 
torus is not fed by an external mass reservoir as in case 
of a collapsar, and the timescale of BH accretion is de- 
termined by viscous transport in the torus, not by the 
collapse timescale of a massive star. The torus is very 
compact. Typically it contains a few hundredth up to 
some 0.1 solar masses and extends to 15-20 Schwarzschild 
radii in the equatorial plane and about half that distance 
in the vertical direction. Its density and accretion rate 
are therefore high so that the torus can become partly 
opaque to neutrinos. With neutrino luminosities in excess 
of lO^'^ergs"^ z^i'-annihilation can deposit energy in the 
close vicinity of the black hole at rates between several 
10^°ergs~^ and more t han lO^^ergs"^ l|Ruffert fc Jankal 
ll99flH.Tanka et alJll999^ . 

The early phase of the propagation of jets 
in coUapsars has been studied by nonrelativi stic 
lMa,cFadven fc Vyoosle^ll 999HMa,cFa,c|ven et aljliool as 



well as spe cial jZ hang fc MacFadvcn 2003") or general 
relativistic l)Alov et al...2000; Alov fc Marti 2002^1 hydro- 
dynamic simulations in two dimensions (i.e., assuming 
axial symmetry). Also flrst axisymmetric magnetohydro- 
dynamic simulations have been done for the formation 
of polar outflows from the ce ntral non-rotating BH in a 
collapsar in both Newtonia n llProga et a l. "200?; using a 
pseudo-Newtonian potential. IPacyiiski fc Wiita 1980.) and 
general relativity ijMizuno et al.l l2004t assuming a fixed 
Schwarzshild metric). Numerical modeling of compact ob- 
ject mergers with some realism of the microphysics input 
(equation of state, neutrino proce sses) has mad e progress 
in case of NS -t -NS b inaries llRuffert et a l.' 'l996a'l^; 
Ruffert fc Jankal bO Olt lR,osswog fc Liebendorfer .2003l 
R.osswQg et al.ll2003l) a i id NS+BH systems iTanka et alJ 
1999t I.Tanka fc Rufi-ert' "2002; ' Rosswog et a1.ll2nn4^ . but 
detailed hydrodynamic simulations of the formation of 
the relativistic outfiow from the accreting post-merging 
BH have not been undertaken prior to the present work. 

In this paper we present results of the first general 
relativistic (a fix Schwarzschild metric is assumed), ax- 
isymmetric models for the acceleration and expansion of 
ultrarelativistic jets and winds, which originate from the 
pair-plasma-photon fireball that is produced by the de- 
position of thermal energy above the poles of a stellar- 
mass BH that accretes matter from a thick accretion 
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torus as obtained in NS+NS or NS+BH m erger sim- 
ulations l)Ruffert k Jankal Il999l: I Janka et all 11999). We 

use the high-reso lution shock-capturing code GENESIS 
llAlov et alJll999t) to integrate the general relativistic hy- 
drodynamic equations in "2.5 dimensions" (i.e., in spheri- 
cal (r, 9) coordinates assuming that there is no azimuthal 
{(/)) dependence of physical quantities although the 0- 
components of vectors can have non-zero values). 

Our simulations attempt to answer questions about 
the coUimation mechanism of the polar outflow, the open- 
ing angle of the ultrarelativistic ejecta, the maximum and 
mean Lorentz factors that can be attained, the inter- 
nal structure of the outflow, the duration of a possible 
GRB event, and the isotropic equivalent energy which an 
observer would infer by assuming the source to expand 
isotropically. 

These questions may at most be guessed but they can- 
not be reliably answered on grounds of merger models 
and a consideration of their energy release by neutrino 
emission and the subsequent conversio n of some of this^ 
energy by t/P-ann ihilation to e"'"e~- pairs (Ruffcrt & Janka 
1999Hjanka et al...l999: iRosswog fc Ramirez-Ruiz .2002 : 
Rosswog et al.ll2003|) . We stress here that the fraction of 
the latter energy which ends up in ultrarelativistic out- 
flow with low baryon loading cannot be deduced or "ex- 
trapolated" from known spatial distributions of mass and 
i^P-energy deposition at any instant of a merger or accre- 
tion simulation without self-consistently time-dependent 
hydrodynamic modeling. The relativistic outflow develops 
in a complex hydrodynamic interaction with the accretion 
torus, cleaning its own axial funnel such that later energy 
deposition encounters a much reduced baryon pollution. 

Answers to the questions mentioned above are of cru- 
cial importance for discussing compact binary mergers as 
potential sources of observed short GRBs and the obser- 
vational properties of those. Current wo rk on these is- 
sues ("e.g-. teosswog fc R,amirez-R,uijl2003j) is much hand- 
icapped by the lack of quantitative information from de- 
tailed numerical simulations. Conclusions in the latter pa- 
per, for example, depend sensitively on the coUimation 
of the ultrarelativistic ou tflow. In this context work by 
iLevinson fc Eichlcr" (200rf) is referred to, who suggest a 
hydrodynamic coUimation of the relativistic fireball by a 
surrounding baryonic wind emanating from the accretion 
torus and calculate the corresponding opening angle by 
analytic solutions. The discussion by these authors, how- 
ever, is based on a number of simplifications and assump- 
tions and therefore awaits verification by more complete 
models. Our simulations, for example, reveal that the ver- 
tical extension of the thick accretion torus cannot be ig- 
nored but plays a decisive role with respect to outfiow col- 
limation, baryon entrainment and structure development 
by Kelvin-Helmholtz instabilities at the jet-torus inter- 
face. 

Our paper is organized as follows. In Sect. 2 we shall 
describe in some more detail our model setup, the initial 
conditions, the numerical methods, and the parameters 
and assumptions which are varied systematically in our 



studies. Section 3 will provide a discussion of our results 
for two types of model sequences with different assump- 
tions about the conditions for the deposition of thermal 
energy in the vicinity of the BH-torus system. We distin- 
guish "type-A" models (in Sect. 3.1) with fairly high den- 
sity in the surroundings of the torus, and "type-B" models 
with low environmental density (Sect. 3.2). In Sect. ^ we 
shall also consider the evolution of type-B models for a 
period of 0.5 seconds after the onset of the energy re- 
lease. In some of these models the energy deposition was 
taken constant up to 0.1 seconds and then switched off. 
In other models the energy deposition was modeled more 
"burst-like" , reaching maximum values for a period of only 
20 milliseconds and then decaying w ith time like a t^'^^'^ 
power law (cf.. lSetiawan et ani2004l) or modulated during 
the whole evolution (including the decay phase) with a 
certain time variability of the energy released per unit of 
time. A summary of our main results and conclusions will 
follow in Sect. 

2. The Model 

Short GRBs (or, at least, a subclass of these events) might 
be produced by the systems resulting from mergers of pairs 
of compact objects (either two neutron stars or a neu- 
tron star and a black hole). The situation that can arise 
from the mentioned merging processes consists of a cen- 
tral black hole of a few solar masses girded by a thick 
accretion torus with mass between 0.05 arid « O.3M0 
(see, e.g..'R, uffert fc Jankalll999t Ijanka et al.l ll999t) . Once 
the thick disk has formed, neutrino-antineutrino annihi- 
lation may release up to ~ 10^^ ergs above the poles of 
the black hole in a region that contains less than 1O~^M0 
of baryonic matter. This may lead to the acceleration of 
this matter to ultrarelativistic speeds, possibly account- 
ing for a GRB with observed properties. If the dura- 
tion of the event is related to the lifetime of the system 
llSari fc Pira nll997^ this kind of events can only belong to 
the class of short GRBs because the expected timescale on 
which the black hole engulfs the disk is at most of a few 
hundred milliseconds jRuffert fc Jankalll999l: Ijanka et alJ 
llQQfltlLee fc R,amirez-R,uijl2002HSetiawan et al.ll2004 . 

In order to test the viability of post-neutron star 
merger systems as likely progenitors of short GRBs we 
have used general relativistic, hydrodynamic simulations 
for studying the formation, acceleration and early propa- 
gation of the relativistic outflow that is driven by the de- 
position of thermal energy near the accretion torus around 
a stellar-mass black hole. In such a black hole-torus con- 
figuration a jet is expected to be launched by any process 
which gives rise to a local deposition of energy and/or mo- 
mentum, e.g., ;/i/-annihilation, or magneto- hydrodynamic 
processes. We mimic this process by releasing thermal en- 
ergy in a prescribed cone around the rotational axis of the 
system. 

The initial model is set up to reproduce the ex- 
pected state of a post-neutron star merger accretion torus 
around a central black hole with a mass in the range 
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Fig. 1. Logarithm of the rest-mass density of the initial model of type- A simulations showing the structure of the 
relaxed accretion torus. The scale of the plot is limited in order to be able to distinguish properly the toroidal 
structure. In the upper left corner the maximum and minimum values of the displayed variable are reported. 



2.44 Mq — 3Mq. Thus, the initial model is not the result 
of a merger simulation of two compact objects (at least 
one of them being a neutron star) but instead, we use two 
different ad hoc procedures to construct it. In the first one, 
a toroidal distribution of matter and angular momentum 
is placed in our computational domain, embedded by a 
larger, essentially spherically symmetric dilute gas cloud. 
S uch a configuratio i i is su ggest ed by the previous results 
of lRuffert &: Jankal l)l999() and I Janka et all l)l999|) where 
the gas cloud is produced by mass shed off the hot massive 
merger remnant. Of course, a more refined setup would re- 
quire to take into account asphericities due to the angular 
momentum in the system. For the time being, however, 
we constrain ourselves to spherically symmetric external 
gas clouds and avoid further refinements. This initial con- 
figuration is not in equilibrium and we need to evolve it 
for some tenths of a second in order to obtain a relaxed 
initial state. The relaxation of the model leads to the shed- 
ding of mass from the torus that expands into the envi- 
ronment. Hence, the relaxed model has a relatively high 
environmental density of ~ 10^ g/cm~^ (Fig- Oil between 
r ~ 10'' cm and r ~ 10^ cm, declining steeply down to 
~ 2.5 • 10^'^g/cm^'^ at r ~ 3 • 10^ cm. We will refer to 
this kind of initial models as type- A models (Fig^). The 
torus and black hole masses of type-A models are 0.17 Mq 
and 3 Mq after the relaxation time, respectively. The inner 
grid boundary is located in type-A models at 4 gravita- 
tional radii {Rg = GM/c^; G, M and c being the gravita- 
tional constant, the mass of the black hole and the speed 



of hght in vacuum, respectively), i.e., at « 1.78 • 10^ cm 
from the center. 

The production of type-B i nitial models u ses a pre- 
scription very close to that of iFont fc Daign3 (|2.002) in 
order to build an equilibrium torus around a Schwarzschild 
black hole. Differences with respect to iFont fc Daignj 
l|2002l) arise from the fact that we do not use a barotropic 
equation of state (EoS) and, therefore, it is not possible 
to obtain from the specific enthalpy distribution (which 
is the primary thermodynamical variable that is obtained 
from the Euler equation) unambiguously all other ther- 
modynamic variables in the equilibrium torus. Thus, the 
adiabatic index (7) at every point has to be guessed, i.e., 
the initial model is not in perfect equilibrium (although 
the initial state is much closer to equilibrium than in the 
previous case), and it also requires some transient period 
until the initial model evolves to a relaxed configuration. 
However, as the deviation from equilibrium is small we 
can start releasing energy right from the beginning of the 
simulation because whatever the changes due to the re- 
laxation are, they are relatively small and they propagate 
radially at a speed smaller than any ultrarelativistic out- 
flow. On the other hand, a neutron star merger remnant 
originates from a violent event and, therefore, both the 
emergent accretion torus and its environment must natu- 
rally be expected to exhibit some variability in time. Thus, 
the initial models of type-B, although not fully relaxed, are 
not in obvious conflict with the situation that is expected 
from merger simulations (FiglSJ. The type-B models have 
a much smaller environmental density than those of type- 
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Fig. 2. Logarithm of the rest-mass density of the initial model of type-B simulations showing the structm'e of the 
accretion torus. The scale of the plot is limited in order to be able to distinguish properly the toroidal structure. In 
the upper left corner the maximum and minimum values of the displayed variable are reported. Note that in contrast 
to Fig. n the torus is much denser and more compact. 



A for radii larger than « 10'' cm and their density gradient 
(in most of the computational domain) is also steeper than 
in the other case {p ~ r"^ ''; Fig.|2I). Models of type-B have 
a torus and a black hole mass of 0.13 Mq and 2.44 M©, 
respectively. The inner grid radial boundary is located in 
type-B models at IRg, i.e., at « 7.57 ■ 10^ cm from the 
center. 

2.1. Numerical details 

In the radial direction the region of energy deposition ex- 
tends from the inner to outer grid boundaries. In the an- 
gular direction, the half-opening angle of the deposition 
cone {0o) around the system axis spans a range which 
was varied between 30° an d 75°. Fro m the results of 
[ Ruff ert fc Janka (1999(1 and I.Tanka elal . (.1999) we have 
inferred a power-law decline for the energy deposition rate 
per unit of volume by i/iz-annihilation in the direction per- 
pendicular to the equatorial plane of the accretion torus. 
We have chosen an explicit form of the energy deposition 
law in terms of observer frame coordinates as 

q = qo{zo/z)", (1) 

where z is the distance along the rotation axis, n is the 
power-law index, zq = 2Rg is the minimum vertical dis- 
tance above the equator where energy is deposited, and 
qo is the normalization factor that we use to set the to- 
tal energy deposition rate E as measured by a static ob- 
server at infinity. In this paper we have considered only 



the case n = 5, because this value complies best with 
the data in Ijanka et al.1 l|l999|) . We assume energy de- 
position rates which vary in the range from lO^^ergs"^ 
to 5 ■ lO^^ergs"^ (for the summed values of both hemi- 
spheres) and which are constant in time. We have also in- 
cluded two cases where the energy release is chosen to be 
time dependent (the explicit time dependence of models 
B07 and BOS is provided in Eas.irUlandlTTl respectively). 
The considered values roughly bracket the expected en- 
ergy deposition rates from z/P-annihilation in the context 
of post-neutron star merger accretion tori. 

The simulations were performed with the multi- 
dimensional relativistic hydrodynamic code GENESIS 
l|Alov et a].lll999h using 2D spherical coordinates (r, 6) . 
We assume equatorial symmetry and we cover the 9- 
direction with 200 uniform zones between 0° and 90°. 
In the r-direction the computational grid consists of 400 
or 500 zones. These are spaced logarithmically between 
the inner boundary located at 4Rg (type-A models) or 
2Rg (type-B models) and an outermost radius of -Rmax = 
3 • 10^ cm (type-A models) or i?max = 2 • 10^° cm (type- 
B models). Even for this moderate resolution, the typical 
number of time steps per model in order to reach a time 
of 0.5 s is « 10^-'^. 

The first set of models (type-A), i.e., those having the 
smaller i?max, is employed to study only the generation 
phase of a relativistic fireball by following the evolution of 
the resulting outflow in a medium of relatively high den- 
sity up to 100 ms. The second set of models, i.e., those hav- 
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Fig. 3. Density vs radius along the symmetry axis of the 
two types of initial models. The dashed-dotted (dashed) 
line corresponds to initial models of type- A (type-B). The 
solid straight lines display power-law fits to both density 
distributions. 



ing the larger i?max, start from non-relaxed type-B initial 
models and, as the computational domain extends up to a 
larger radius, allows us to study the evolution of the fire- 
balls in relatively low-density environments until ^ 0.5 s. 
We have also extended the domain of one of the most 
energetic type-A models in order to follow the evolution 
of fireballs in relatively high-density environments until 
- 0.4s. 

The space-time around the non-rotating black hole is 
described by the steady spherical Schwarzschild metric. 
Effects on the dynamics due to the self-gravity of the ac- 
cretion torus or the external gas are neglected, i.e., we 
consider only the gravitational potential of the black hole. 
The values of the energy deposition rate that we will use 
throughout the paper are measured by an observer at in- 
finity. Therefore, the energy deposition rate per unit of 
coordinate volume q measured by such an observer has 
to be blue-shifted in order to determine the local value 
of the energy deposition rate in the close vicinity of the 
BH. This value has to be further transformed to the local 
comoving fluid frame because it is numerically included in 
the algorithm as a source of the comoving energy density. 
The amount of energy released per unit of coordinate time 
and proper volume in the comoving frame is 



(2) 



" (l-2i?g/r)i/2(^(r2-l) + l) ■ 

where T is the bulk Lorentz factor of the moving fluid 



The EoS assumes nuclei to be disintegrated to 
free, non-relativistic nucleons, treated as a mixture of 
Boltzmann gases, and includes the contribution from ra- 
diation, and an approximate corre ction due to e+e~-pairs 
as described in lWitti et alJ l|l994l) . Complete ionization is 
assumed, and the effects due to possible degeneracy are 
neglected. 
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Fig. 4. Logarithm of the rest-mass density (top panel), of 
the pressure (middle panel) and of the fluid Lorentz factor 
(bottom panel) vs radius along the symmetry axis for the 
models AOl (solid lines), A02 (dashed lines), A05(dashed- 
dotted lines) and A09 (long-dashed lines) after 100 ms. 
The displayed time is larger than the transition time in 
models AOl and A02 and, therefore, the recoUimation of 
the initial ultrarelativistic wind solution yields relativistic 
jets in these two cases. 



3. Results 

Given the fact that our models do not include a fully con- 
sistent treatment of all relevant physics and our initial con- 
ditions are set up ad hoc instead of resulting from merger 
simulations of compact binary stars, we have performed a 
series of simulations trying to cover a physically meaning- 
ful parameter space. Nonetheless, we point out here that, 
even if we had done consistent simulations, the variety of 
resulting initial configurations could be very large because 
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the properties of the forming black hole-accretion torus 
system and its environment are sensitive to a number of 
parameters which are incompletely known or variable (the 
masses and spins of the binary objects, the neutron star 
equation of state, general relativistic effects in the binary 
merging, etc.). 

According to iRuffert fc Jankal l)l999|) . the accretion 
torus might be surrounded by a halo of relatively high 
density (a by-product of the merger) that could extend 
up to a few hundred kilometers from the center of the 
merger. On the other hand, merger events should mostly 
occur in the interstellar medium and, therefore, a merger 
and its halo should be surrounded by a very low-density 
environment. The size and the exact density profile of 
the halo can strongly depend on, e.g., the physical pa- 
rameters of the merging progenitors, the details of the 
merger dynamics, the time delay of the black hole for- 
mation, and the neutrino-driven baryonic outflow from 
a transient, massive, hot post-merging neutron star. A 
high density halo would have to be expected if, for ex- 
ample, the collapse to a black hole were delayed due to 
the effects of yerv rapid (differ e ntial) rotation or visc ous 
heating fe.g.. IPuez et aTll2nn4 iMorrison et a].ll2nn4|>. in 
which case the hot neutron star w ould radiate neutri- 
nos and a neutrino-driven wind fe.g.- lDuncan et aDll986|) 
would lead to a dense, expanding baryonic cloud around 
the merger site. It may also be possible to find situations 
where the accretion torus is surrounded by a thin, dilute 
halo, in particular, if the BH forms during the merger 
or within a few dynamical timescales afterwards (e.g., 
Shibata k Unill200d: IShibata et al]l2003t lOechsUn et alJ 
20041) . In regard of these considerations we have performed 
two series of simulations. The first one (models AOl to 
A09 of Table ^1 uses an initial model of type-A and cor- 
responds to a case where the density of the halo is high 
and the density dechne with distance is moderate (Fig.|3I). 
The second set of models (Bl to B8) start from an initial 
model of type-B and, therefore, they correspond to the 
low-density halo class. 

3.1. The case of a high density halo 

For models of type-A, the parameter study comprises 
three different aspects. The first one is the study of the 
effects associated with a variation of the total energy de- 
position rate from E — 10'*^ ergs~^ to i? = 5 ■ 10^^ ergs~^, 
fixing the value of the opening half-angle of the depo- 
sition region at 9q = 30° and using a power law index 
of the energy deposition profile of n = 5. The models 
computed for this purpose are AOl, A02, A05 and A09 
(Table nj. The second aspect is the dependence of the re- 
sults on 00 while keeping the rest of the parameters fixed 
{E = 2-10^° ergs~^, n — b). The models used in this study 
are A02, A03 and A04 (Table HI). Finally, we have made 
two experiments where the energy released per unit of vol- 
ume and time (the deposition rate density, hereafter) was 



kept constant but both E and were changed (models 
AOS and AOS). 

As a general conclusion, we find that, for energy depo- 
sition rates larger than a certain threshold, all the models 
lead to a relativistic or ultrarelativistic outflow (i.e., a fire- 
ball) . This threshold is caused by the need to overcome the 
ram pressure Pram that the external medium exerts on the 
nascent fireball close to its initiation site. If the amount 
of energy per unit of time and per unit of volume that is 
pumped into the region extending up to the scale hight 
of the torus {H ~ lO*" cm; note that about 97% of the 
energy is released in this region due to the large value of 
n in Eq. ^) in a time interval of about the free-fall time 
of fluid elements at a distance of the order of H , is less 
than Pram, then the fireball is unable to expand, and is 
swallowed by the BH. The precise value of this thresh- 
old is model dependent. The threshold is higher for halos 
of higher density than for dilute ones. The shape of the 
fireball and its dynamical features depend very strongly 
on the type of initial model that is considered as will be 
discussed below. 

3.1.1. Dependence on the total energy deposition rate 

As we have already pointed out, there is a threshold en- 
ergy deposition rate below which a relativistic fireball 
is not initiated. For models of type-A, the threshold is 
Ethv ^ 1049ergs-i(for Bq = 30°), while for models of 
type-B Ethv < 10*^ ergs"i(for 6*0 = 45°). The difference 
between both thresholds comes from the fact that the 
higher halo density of models of type-A leads to a higher 
ram pressure than in models of type-B. 

In models of type-A, when E is above i^thr, the result- 
ing outflow consists of either (i) a stratified (both in radial 
and polar directions) wind of approximately conical shape 
or (ii) a relativistic jet. Let us note that, by definition, the 
difference between a wind and a jet is that in the latter 
the lateral boundaries are causally connected while they 
are not so in case of a wind. This distinction is relevant 
when the stability properties of both kind of solutions are 
considered^. 

A wind results for high energy deposition rates per 
solid angle (models A05 and A09, Figs.ElEl)- In the polar 
direction, the wind is confined in a cone the angle of which 
is model dependent (see below) . Its structure consists of a 
sheared wind with a central, low density, unshocked, ultra- 
relativistic part (near the symmetry axis) surrounded by 
a slow (mildly relativistic), denser, shocked layer. Surfaces 
of constant physical variables are almost spherical in the 
ultrarelativistic part of the wind. In the radial direction, 
the wind extends radially from a fiducial point (located 
at ~ 3 • 10^ cm) to a terminal shock where it interacts 
with the external medium. Due to the larger density of 

^ For example, a jet can be destroyed by Kelvin Helmholtz 
instabilities that originate at its surface, while a wind only 
develops a turbulent layer which does not affect the whole so- 
lution. 
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Table 1. The first three columns show the model name and the two main parameters which characterize the energetics 
of the models, namely the energy deposition rate E (which represents the summed value for both hemispheres) and 
the half-opening angle of the cone over which the energy is released, respectively. For the models B07 and BOS entry E 
gives the peak values of the energy deposition rate {Eq, Eg. Columns four, five and six give some kinematic and 
dynamic properties of the outflow: the velocity {vp) of the propagation of the fireball in the radial direction (i.e., the 
velocity of the leading, radial edge of the outflow), the maximum Lorentz factor (Fmax) attained in the outflow and, 
the wind (or jet) half-opening angle (6*^1) of the outflow after 10 ms or 100 ms, respectively. The last column gives the 
mass of the outflow (M/) after 100ms. For models A04 and A06, there are data only up to 10ms. The half-opening 
angle of the outflow is computed as the maximum ^-coordinate of all computational cells in where the fluid Lorentz 
factor exceeds a value of 10 and the radial velocity is positive. The mass of the fireball is computed by adding up the 
mass in all computational cells that match the same criterion as the one used to compute 9^ ■ The choice of a relatively 
small value of the Lorentz factor (F = 10) in our criterion to separate the fireball matter from the external medium 
is motivated by the still ongoing acceleration of the fluid in the fireball that takes place at the expense of the large 
internal energy of the outflow. 
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the external medium, two shocks form: a reverse and a 
bow shock. Through the reverse shock the wind deceler- 
ates from ultrarelativistic velocities to mildly relativistic 
ones (between 0.6c and 0.97c, see Vp in TableCJ. The prop- 
agation speed of the bow shock is smaller than that of the 
fluid in the wind (Fig. 0)). This explains why the radial 
expansion of the wind is only moderately relativistic. 

In a conically expanding relativistic wind the physical 
variables change with the distance z along the axis of ro- 
tational symmetry accor ding to simple power laws (e.g., 
iLevinson fc Eichleill200o|) . For example, the pressure, the 
density and the Lorentz factor of a fluid whose adiabatic 
index is 7 = 4/3 follow the relations: 



p/Po^{z/zo) ^, pIpo^{z/zo) ^, r/To^z/zo, (3) 



where pq, po and Fq arc the values of the pressure, of the 
density and of the Lorentz factor at the fiducial point of 
the wind located at a distance zq. In our models (Fig. 0)), 
we obtain for the unshocked region of the outflow (the 



shocked one does not follow a power law) 

p/po = {z/zo)-"^-'', (4) 
p/po = (^/^o)~'■'^ 

r/To^{z/zor''. 

The discrepancy in the power-law indices arises from the 
fact that the value of 7 is slightly smaller than 4/3 in the 
wind of our models due to the presence of e''"e~-"pairs in 
the fluid. Since for type-A models the pressure and the 
density in the wind decrease faster than in the external 
medium up to 2 10^ cm (Fig. ^ the outflow turns into 
a jet, because the wind pressure eventually becomes equal 
to that of the external medium. At that point, a recol- 
limation shock forms starting from the wind surface and 
propagating towards the axis. Since the fiducial values of 
the pressure increase with an increasing E, a higher E 
reduces the possibility of matching the pressure of the ex- 
ternal medium. Therefore, the structure of the wind tends 
to remain conical when the energy deposition rate is in- 
creased. If the recollimation shock does not form before 
the fireball reaches the part of the environment where a 
steep density decline sets in ( ^ 10^ cm; Fig. OJ, it might 
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Fig. 5. Logarithm of the rest-mass density for models AOl (panel a), A02 (panel b), A05 (panel c) and A09 (panel d) 
after 100 ms of energy deposition. The color scale is limited in order to enhance the details of the outflow region. The 
maximum and minimum density values (in logarithmic scale) of each model are given below the panels. 



not form at all (because for radial distances larger than 
~ 10^ cm the rate of decline of the pressure in the external 
gas exceeds that in the ultrarelativistic wind). 

When the energy deposition rate is only slightly 
(model AOl) or moderate (model A02) above the thresh- 
old, the conical wind converts to a relativistic jet due to 
the recoUimation shock (Fig. O. The outflow then pos- 
sesses all the morphological elements known for light rel- 
ativistic jets, i.e., a knotty, unsteady beam that ends in a 
Mach disc surrounded by an over-pressured cocoon. The 
outflow begins as a transient, expanding, over pressured, 
conical wind. But after a transition time ( Attrans ~ 10 ms 
for model AOl and Attrans ~ 30 ms for model A02) a rec- 
oUimation shock forms due to the pressure confinement of 
the external medium reducing the opening angle to val- 
ues less than 1°. The formation of the recoUimation shock 



has two main consequences. First, it separates the outflow 
in two parts: an unshocked inner flow (close to the depo- 
sition region; the smooth region in Fig. and an out- 
her shocked flow where the pressure and the density have 
substantially larger values (Fig.0J. Second, the recoUima- 
tion shock effectively stops the linear increase (Eq. O of 
the beam Lorentz factor observed in the more energetic 
models (A05 and A09; see Fig.0J) even at large distances. 
Hence, the propagation speed of the head of the outflow 
is about 35% smaller in models AOl and A02 than that in 
models with higher energy deposition rates (Table ^| af- 
ter the initial transient phase. As the fiducial pressure of 
the fireball grows with E, while the pressure decrease with 
distance is approximately the same for the wind outfiow 
in all models (Eq. , the distance where the pressure of 
the wind and the external medium match increases with 
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Fig. 6. Logarithm of the Lorentz factor for models AOl (panel a), A02 (panel b), A05 (panel c) and A09 (panel d) 
after 100 ms of energy deposition. The color scale is limited in order to enhance the details of the outflow region. The 
maximum and minimum density values (in logarithmic scale) for each model are given below the panels. 



growing E (the same trend holds for Aitrans)- Whether 
pressure matching occurs or not depends on the value of 
the fiducial pressure and on the pressure gradient in the 
halo. Within 100 ms this happens in none of the more en- 
ergetic models (A05 and A09). 

Increasing the energy deposition rate in models of 
type-A leads to a progressive widening of the resulting fire- 
ball (Fig. O . For models that produce an ultrarelativistic 
wind (but not a relativistic jet), Tabled shows that the 
half-opening angle of the wind is smaller than the half- 
opening angle 6*0 of the energy deposition cone (in case of 
a jet the half-opening angle is set by the interaction with 
the external medium and the cocoon dynamics). Thus, the 
opening angle is not determined by our ad-hoc choice of 
^0- Instead, the opening angle is defined by the angular 
coordinate where the pressure of the fireball equals that 



of the torus. This angle depends on two factors: (i) the 
torus pressure at a given radial distance, increases as one 
moves from the axis to the equator and, (ii) the pressure 
of the fireball, at a given radial distance, is an increas- 
ing function of the energy deposition rate. Hence, as we 
increase E the location where the pressure of the fireball 
matches that of the torus shifts to larger angles, i.e., the 
opening angle of the outflow grows. 

Within an outflow crossing time of the torus in the ax- 
ial direction ( ^ 1 ms), the fluid accelerates up to Lorentz 
factors of w 10 (Fig. and, if it does not form a jet, the 
fluid becomes ultrarelativistic, i.e., any motion except in 
the radial direction is prevented. Hence, a sideways ex- 
pansion is drastically reduced, and the initially imprinted 
opening angle remains unchanged. After 10 ms the max- 
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imum Lorentz factor rises to ~ 100 in the case of winds 
while it is only ~ 18 for jets. 

The amount of mass in the fireball (Table defined 
as the sum of the masses of each computational cell in the 
outflow that moves out radially with a F > 10, increases 
almost linearly with E during the period in which energy 
is released in the system. The mass is mainly concentrated 
in the cocoon that surrounds either the outflowing jet or 
wind. In the case of winds, only a small fraction of the 
mass moves at Lorentz factors larger than 100. 

3.1.2. Dependence on the half-opening angle of the 
energy deposition cone 

We have performed a series of runs where we compared 
for fixed energy deposition rates the influence of an in- 
creasing angular width of the region of energy release, Oq 
from 30° to 75°. The models considered in this study are 
divided into two groups. The first group consists of A02, 
A03 and A04 (Table [H Fig. d) which have a moderate 
energy deposition rate of E = 2 • 10^° erg/s (group M 
hereafter). The second group possesses a higher deposi- 
tion power of 10^^ erg/s, and consists of models AOS, A06 
and A07 (Table H) Fig.il group H hereafter). 

For a flxed, moderate energy deposition rate the open- 
ing angle of the outflow is initially quite insensitive to the 
increase of the deposition angle Oq: for the M-group, while 
00 grows from 30° to 75°, O^j only decreases from 11.3° 
to 8.5° (we define as the angle of the outflow hav- 
ing F > 10) in 10 ms. After longer evolutionary times all 
the models of group M tend to form relativistic jets with 
opening angles that depend on the interaction between 
the beam, the cocoon and the external medium, but not 
on 9q. For models of group H, depending on the value of 
^0, we obtain outflows with opening angles after 100ms 
of evolution that exhibit a rather small range of variation 
between 15° (A05; 0o = 30°) and 8.13° (A07; 0o = 75°). 
Hence, the trend for group H is a slight non-monotonic 
decrease of the outflow opening angle as we deposit the 
energy in a wider region. 

These results confirm, on the one hand, that the open- 
ing angle of the deposition region is not directly related 
to the final opening angle of the outflow and, on the other 
hand, that the mechanism of coUimation of this series of 
models is significantly different from the mechanism dis- 
cussed bv lLevinson fc Eichleil l|200(]l) . These authors pro- 
pose that a baryon-rich wind emerging from the torus is 
able to coUimate the baryon poor jet moving along the 
symmetry axis of the system. For type-A models the role 
of such a baryon-rich wind is played by an almost static 
cloud that is formed during the initial relaxation phase of 
the torus. This cloud has a rather low velocity (~ 0.027c) 
and a density of ~ lO^g/cm^^ at a distance of 10^ cm. 
The density does n ot fall off as (as in c ase of the 
baryon-rich wind of iLevinson fc Eichleil l200o|) . but with 
a much smaller power. Therefore, the main effect of the 
environment is, in our case, to provide sufficient inertial 



confinement for preventing excessive sideways expansion 
of the axial ultrarelativistic flow. 

The insensitivity of 0^, to ^o, arises because the open- 
ing angle of the relativistic outflow is mainly determined 
by the inclination angle of the torus walls around the rota- 
tion axis. As discussed in ? 13. 1.11 the torus does not have 
a sharp surface that separates it from the halo. Instead, a 
gradient of density and pressure connects it with the sur- 
rounding medium. We have argued that the region where 
the pressure of the torus and of the emerging fireball 
match each other determines the initial 0^. In the con- 
sidered set of models, the energy deposition rate per unit 
of volume decreases with growing ^o- This fact implies a 
smaller pressure (at any given radial distance). Thus, pres- 
sure equilibrium between the fireball and the torus occurs 
at a smaller polar angle, which in turn reduces the value 
of 0w as can be seen in Figs. and |H1 

Another consequence of increasing 0q is the transition 
from a sheared, ultrarelativistic wind to a knotty, moder- 
ately relativistic jet. For models of group M and 6*0 > 45°, 
this transition takes place very soon after the birth of the 
fireball (A^trans ^ 10 ms; Fig. d) at small radial distances 
(models A03 and A04). In the case of the largest 0o con- 
sidered (model A04, 0q = 75°), the maximum Lorentz 
factor is only 6.97, and the jet is prone to large ampli- 
tude Kelvin-Helmholtz (KH) instabilities that lead to a 
wrinkled surface and an inhomogeneous beam. Increasing 
the energy deposition rate (models of group H), we ob- 
serve the same qualitative behavior, but the wind-to-jet 
transition happens at a larger value of do (« 75°). 

Model A04 produces an inhomogeneous, low- velocity 
jet because the imposed energy deposition rate per unit 
volume for the chosen value of is only slightly above the 
threshold value of E that is necessary to launch an out- 
fiow at all. We have confirmed this by performing another 
simulation with the same energy deposition rate but with 
6*0 — 90°, which did not produce any outflow. For a given 
threshold value for the energy deposition rate per unit of 
volume, the minimum total energy deposition rate to pro- 
duce an outflow increases with larger 00 . The formation 
of jets or winds depends on the energy released per unit 
of time and of volume. The larger the energy release, the 
larger is the chance for wind formation. 

3.1.3. Dependence on the deposition volume Vjcp and 
on the total energy deposition rate E for 
-E/Vdep=constant. 

In this section we will study the dependence on the depo- 
sition volume if the rate per unit volume at which energy 
is released in the system is fixed. The two models that we 
consider here are AOS (having a large deposition rate and 
a large deposition angle) and A08 (with smaller deposition 
rate and deposition angle but the same energy released per 
unit of time and volume as A03). Both models develop a 
mildly relativistic jet and display very similar properties 
at the beginning of the evolution (Fig.EJ. However, they 
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Fig. 7. Logarithm of the Lorentz factor for models A02 (panel a), AOS (panel b) and A04 (panel c) after 10ms of 
energy deposition. The color scale is limited in order to enhance the details of the outflow region. The maximum and 
minimum Lorentz factor values (in logarithmic scale) of each model are given below the panels. 



later evolution differs. For example, the propagation veloc- 
ity of the jet head is smaller for AOS than for AOS, because 
model AOS accumulates more mass in the cocoon, and a 
part of this mass enters the beam thus leading to a more 
massive jet {Mf°^ = 4.5 • lO^^ g while Mf°^ = 1.4-1022 g). 



respondingly lower, although the rate of energy deposition 
per unit volume is the same as in model AOS. 



This can be understood from the fact that in model 
AOS the cone where energy is deposited has a larger half- 
opening angle 6^ and thus overlaps more with the outer 
layers of the torus than it does in model AOS. The en- 
ergy release in this overlap region causes mass shedding 
from the outer layers of the torus where the densities (for 
SI. 4° < 9 < 45°) are 2 to 4 orders of magnitude higher 
than those in the solid angle corresponding to the deposi- 
tion cone of model AOS. Hence, the outflow in model AOS 
becomes baryon-enriched. The outflow velocities are cor- 



Independent of the differences in mass, velocity and 
maximum Lorentz factors reached in the two models, the 
final opening angle becomes less than 4° in both cases 
(TablcQ, because a series of recoUimation shocks in the 
beam prevents a large sideways expansion of the forming 
jet. Nonetheless, the cocoon of the jet of model AOS is 
twice as thick (w S • 10^ cm when the jet has propagated 
for 100 ms) as that of model AOS. These differences in 
the cocoon can be understood from the different angular 
extensions of the deposition regions, too. 
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Fig. 8. Logarithm of the Lorentz factor for models A05 (panel a), A06 (panel b) and A07 (panel c) after 10ms of 
energy deposition. The color scale is limited in order to enhance the details of the outflow region. The maximum and 
minimum Lorentz factor values (in logarithmic scale) of each model are given below the panels. 



3.1.4. Evolution after the switch off of the central 
energy source 

We have done simulations where, starting from model A09 
(which has the largest E of all type- A models), the energy 
deposition was switched off after a source activity time 
tsa. = 100 ms. In this case the ultrarelativistic wind so- 
lution was destroyed within less than 2isa, because the 
fireball sweeps up mass when it moves through the high 
density halo in models of type-A fFigs. [Tnilll|l . The swept 
up mass prevents the terminal shock of the ultrarelativis- 
tic wind from moving at ultrarelativistic speeds (it actu- 
ally propagates at a speed of « 0.97c at 100 ms). Therefore 
the ultrarelativistic wind is decelerated when it catches up 
with the reverse shock fFigs.irUb.lllb). In case of type-A 
models that form jets, a similar argument can be used to 
explain the mildly relativistic speeds at which the heads 



of the jets propagate (w 0.67c). Hence, models of type-A 
will, in general, not lead to successful GRBs, because they 
accumulate too much mass in the fireball and decelerate 
well before the outflow becomes optically thin. 

Note that a necessary condition for a fireball to pro- 
duce a successful GRB is that the rest mass carried by 
the fireball multiplied by c^ should be about two orders 
of magnitude smaller than the energy of the fireball by 
the time at which transparency sets in. If the mass of the 
fireball grows too much before the fireball reaches the pho- 
tosphere, the Lorentz factor corresponding to the terminal 
shock of the wind, and as a consequence, the fiuid Lorentz 
factor, will be reduced below the critical limit (namely, 
F « 100) to account for the properties of observed GRBs. 
An estimate of the maximum Lorentz factor the fireball 
could ultimately reach (if all the internal energy could 
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Fig. 9. Logarithm of the rest mass density (right two panels) and of the Lorentz factor (left two panels) for models 
AOS (panels a), and AOS (panels b) after 100 ms of energy deposition. The color scales are limited in order to enhance 
the details of the outflow region. The maximum and minimum values of the respective variable are given below the 
panels. Note the strong pinching of the jet of model AOS at r « 10* cm which is due to the shedding of the most 
external shells of the accretion torus. 



be converted into kinetic energy) is given by the ratio 
r] = E/Mc^, where E is the total energy of the fireball at 
any given time after the energy deposition is switched off. 
M is the mass of the fireball computed according to the 
definition given in the caption of Table In model A09 rj 
decreases with time once the energy deposition is switched 
off, because of the increasing mass of the fireball. Thus, 
this model will not meet the criterion rj ss 100 which is 
required for producing an observable GRB. 

Although type-A models are not expected to produce 
GRBs, we can estimate the kind of observational signature 
that this kind of events might generate. We assume that 
after some hundred milliseconds of evolution the bipolar 
outfiow can be roughly approximated by two spherical, 
homogeneous clouds of hot radiation-dominated plasma, 
each with a mass M ~ 1O^^M0, radius Rq ^ 10^ cm, 
temperature Tq ^ 1-5 • 10^ K and internal energy Uq ^ 
2 • 10^*^ erg (Figs.EJiEt)- Assuming that the subsequent 



evolution of such a cloud is adiabatic, its temperature 
and internal energy will decrease as the cloud expands. 
Transparency is reached when the radius of the plasma 
cloud becomes 

1/2 / ^ X 1/2 



i?t « 3 • 10^^ cm ( — 



(5) 



_10-5Mq^ 

on a time scale (assuming that the cloud expands with 
about the light speed) 

1/2 / X 1/2 



tt « Rt/c « 10" s — 



/ \ 10-5 Mq 



, (6) 



where k « Kg ~ 0.2cm^g-^ is the mean opacity (with 
Ke being the opacity caused by electron scattering for an 
electron-to-baryon ratio of 0.5). 

When transparency sets in, both the temperature and 
the internal energy of the cloud have decreased to the 
following values 
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(a) A09: 1 = 0. 130 s (b) ACQ: t - 0.200 s (c) A09: t = 0.383 s 
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Fig. 10. Snapshots of the radial velocity of model A09 after the shutdown of the energy release (the snapshot times 
are provided above each panel). The color scale is limited in order to enhance the details of the outflow region. The 
maximum and minimum radial velocity values of each model are given below the panels. Although the outflow has 
reached Lorentz factors above 700 after 100 ms of evolution, and continues accelerating (panel a) until the rear end of 
the outflow catches up with the reverse shock (panel b), after 500 ms only a subrelativistic outflow has survived and 
propagates with radial velocities of ^ 1.25 



Tt « 5 • 10^ K 



To X/.X-/^ 7. 10- ergs- ^^^^ 



1.5- io9k; v«^ey / M vV ^0 



(9) 



10 ^MqJ V 109 cm/' Therefore, we expect a low-luminosity, soft UV-flash 

to be emitted as a result of the nonrelativistic outflow 
expanding from BH-tori systems in high-density merger 
halos. Due to the small value of L.,„ only galactic events 



Ut « 7 - 10*^ erg cm 



To f \ might be detectable. 



1.5- 109 

M V'^' f Ro V 3.2. The case of low-density halo 



10-5 Mp) / \10^< 



(8) 



In case of a type-B progenitor, our parameter study com- 
prises the same three different aspects as for type-A mod- 
After optically thin conditions are reached, most of the els. We first study the dependence on the energy deposi- 
internal energy of the cloud will be radiated over a time tion rate which is varied from lO''^ ergs^^, to 2-10^° ergs~^ 
scale of the order of U and, therefore, the peak luminosity keeping the value of the half-opening angle of the depo- 
w Ut/tt will be roughly sition cone fixed to ^o = 45° and the power-law index of 
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(a) A09; t = 0.130 s (b) A09: t = 0.200 s (c) A09: t = 0.383 s 
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Fig. 11. Snapshots of the logarithm of the rest mass density of model A09 before the shutdown of the energy release 
(the snapshot times are provided above each panel). The color scale is limited in order to enhance the details of the 
outflow region. The maximum and minimum rest mass density values (in logarithmic scale) of each model are given 
below the panels. 



the vertical decrease of the energy deposition rate set to 
n = 5 (see §|OJ. The corresponding models are BOl, B04 
and BOS (Table The second aspect is the dependence 
of the results on 9q ior E = 2 ■ 10^°ergs~^ and n ~ 5 
(models BOl, B02 and BOS; TableP). Thirdly wc compare 
models where the energy deposition rate per unit volume 
is fixed, but was varied (models B02 and B06). 

We have further computed two models where the en- 
ergy deposition rate was assumed to vary with time (mod- 
els B07 and BOS). These models were set up according to 
the following prescription for the energy deposition rat e 
(guided by the results obtained bv lSetiawan et aLll2004j) 



B07 



Eq 
Eq 

Eq 



t 



0.01 

t 



0.03 



-3/2 



for t < 0.01 s, 

for 0.01s < t < 0.03 s, 

for t > 0.03 s. 



(10) 



where Eq = 2.35 • lO'^^ergs'^ With this value of Eo the 
same amount of energy is released in model B07 until time 
infinity as in model BOl during its on-time of 100 ms. We 
point out that due to the rather steep power-law decrease 



of the energy deposition rate with time in B07, more than 
60% of the total energy is released within the first 100 ms. 
In model BOS, the energy deposition rate follows the same 
average time dependence as in B07 but is sinusoidally 
modulated with a time period tvar ~ 10 ms according to 



Ebos = Eboi ■ (0.5sin(27ri/ivar) + 1) ■ 



(11) 



In all other models of type-B we have initiated the 
energy deposition by raising the energy deposition rate 
linearly from zero to its terminal value within 10 ms. This 
timescale leads to a more modest onset of the hydrody- 
namic evolution. It can be motivated by the fact that in 
neutron star mergers the neutrino luminosities and, there- 
fore, the energy deposition by vD annihilation increase 
within a few milliseconds in response to the viscous shear 
heating of the merger remnant (se e, e.g., Ijanka et alJ 
ll999HRosswog fc LiebendorfeJl2003() . 

In type-A models, a relativistic or ultrarelativistic out- 
flow (i.e., fireball) is obtained when the energy deposition 
rate is larger than a certain threshold value I3.1|l . The 
threshold for type-B models is much smaller than that for 
models of type-A in case the opening angle of the depo- 
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Table 2. Type-B progenitor models. The first column gives the model name, and the second one the total energy Ed 
released in two hemispheres during the simulated time. Columns 3 to 6 display the ratio of the total energy (internal 
plus kinetic) in the region of the outflow that has a Lorentz factor larger than a certain threshold value (given by the 
subscript of the energy in each column) relative to the energy that is released (in one hemisphere) . The seventh column 
provides the fraction of the deposited energy (per hemisphere) that is converted into kinetic energy of the outflow 
with Lorentz factor larger than 100. Columns 8 to 10 show the maximum Lorentz factor Fmax, the half-opening angle 
Ow and the mass of the fireball Mf, respectively. The half-opening angle of the fireball is computed as the maximum 
^-coordinate of all computational cells where the fluid has a Lorentz factor above 10 and positive velocities. The mass 
of the fireball is computed by adding up the mass in all computational cells that match the same criterion as the 
one used to compute 9w The last column displays an estimate of the asymptotic value of the Lorentz factor of the 
outflow computed as Foo — Et>w/ [MfC^), i.e., assuming that all internal energy is eventually converted into kinetic 
energy, and no mass is subsequently swept up from the environment. All the data in the table correspond to a time of 
evolution of 0.5 s. 



Model 


Ed [erg] 


-Br > 100 
Ed 


-E'r>50 
Ed 


-Er>io 
Ed 


Et>2 

Ed 


-E*:,r>ioo 
Ed 


Tniax 


dm 


Mf 


[g] 




BOl 


2 • 10« 


0.29 


0.36 


0.59 


0.70 


4.4 • 10"^ 


859 


24° 


7.4- 


10^* 


1765 


B02 


2 ■ 10*^ 


0.09 


0.15 


0.35 


0.48 


4.8 • 10"^ 


687 


15° 


6.4 ■ 


10^* 


1217 


BOS 


2 . lO-"^ 


0.00 


0.00 


1.2 • 10"* 


0.04 


0.00 


16 


3° 


1.9 ■ 




142 


B04 


10« 


0.17 


0.26 


0.52 


0.65 


5.1 ■ 10"^ 


492 


15° 


3.6 ■ 


^q23 


1601 


BOS 




0.30 


0.37 


0.60 


0.71 


8.3- 10"^ 


979 


25° 


3.6 ■ 


^q25 


1848 


B06 


10« 


0.19 


0.29 


0.54 


0.66 


3.7-10"^ 


717 


18° 


3.4- 


10^* 


1761 


B07 


1.65 ■ 10"^ 


0.14 


0.25 


0.59 


0.72 


3.4 • 10"^ 


839 


25° 


7.6 ■ 


10^* 


1429 


BOS 


1.67 • 10"^ 


0.18 


0.25 


0.49 


0.60 


4.5 • 10"^ 


789 


21° 


5.6- 


10^* 


1607 




7.7 8.6 9.4 
Log 10 R [cm] 



10.3 



Fig. 12. Logarithm of the rest-mass density (top panel; 
in gcm~'^), of the pressure (middle panel) and of the fluid 
Lorentz factor (bottom panel) vs radius along the symme- 
try axis for the models BOl to B06 after 100 ms. 



sition cone is chosen to be the same (for — 30° the 
threshold is well below 10''® erg s~^). This can be under- 
stood by the lower ram pressure that the low-density gas 
exerts on the fireball close to its site of initiation in type-B 
models. 

The shape and the internal structure of the fireballs 
in models of type-B (e.g.. Figs. I15land ll6|l differ substan- 
tially from those of type- A models (e.g., Fig.lHI). Fireballs 
in type-B models that reach ultrarelativistic speeds (all 
the B-cases shown in Tabled except BOS) display a con- 
ical core with an opening angle of « 5° to w 20° which 
spreads into a 26° to 35° wide conical region at larger 
radii. In addition, the full structure is surrounded by an 
approximately spherical shell of matter that expands with 
moderate Lorentz factors (< 10). In this shell the variation 
of the Lorentz factor with polar angle depends strongly on 
the model. This structure is produced by the deposition 
of energy in the surface layers of the torus. Different from 
models of type-A, the unshocked region of the outflow (i.e., 
the region of the outflow between the fiducial point and 
the first recoUimation shock) extends to small radial dis- 
tances (after 100 ms, this region covers ^ 10^ cm; Fig. 1121 
except for models B02 and B03) and there is no simple 
stratification in the angular direction. In models of type- 
B KH instabilities originating in the walls of the torus 
are stronger. They produce variations of Lorentz factor 
fFig. ll3|) and density fFig. ll4|) in the outflow, in particular 
close to the angular edge of the central conical structure. 
In the radial direction the KH instabilities cause a mod- 
ulation of the on-average increasing Lorentz factor, and 
local fluctuations of the density and the pressure. 

The growth of KH modes from the interface between 
two relativistic fluids depends, among other parameters, 
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(a) B01 ; t = 0.005 s (b) B01 : t = 0.010 s (c) B01 : t = 0.050 s 
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Fig. 13. Snapshots of the logarithm of the Lorentz factor of model BOl before the shut off of the energy release (the 
snapshot times are provided above each panel). The color scale is limited in order to enhance the details of the outflow 
region. Note that the distance scale has been adjusted in the two left panels in order to enlarge the details of the 
earlier phases of the evolution. The maximum and minimum Lorentz factor values (in logarithmic scale) of each model 
are given below the panels. 



on the contrast of density (x) and of specific enthalpy 
between the two fluids. Roughly, the larger the density 
contrast is, the faster the growth of the KH modes. Thus, 
the faster growth of KH instabilities close to the accretion 
torus in most of the models of type-B can be attributed to 
the larger density contrast between the torus and the out- 
flow (x = ptorus/poutflow — 10"^'' ) as compared to models 
of type-A (x — 10^^). As a consequence of the growth of 
KH instabilities, mass is entrained in the relativistic out- 
flow which modifies the local speed of the fiuid. The to- 
tal mass entrained is comparable to that found for type-A 
models. However, in contrast to those models, the amount 
of matter piled up in the front part of the outfiow is much 
smaller now (because the halo is less dense). This allows 
a highly relativistic propagation velocity of the fireballs of 
type-B models (Table P). Indeed, the propagation speed 
is so large that models of type-B are viable to produce 
GRBs as will be further discussed in § 01 

The mixing of non-relativistic baryons with the rela- 
tivistic fiuid of the outfiow leads to a change of the adi- 
abatic index in the fireball which occurs mainly along 



the boundaries of the outfiow. In a cone of about 20° 
the outflow is still strongly radiation dominated because 
the average temperatures and densities are 3 • 10'' K and 
10~i2 g cm^-^ (except in case of models B02 and BOS where 
the density is much larger, although they are still radia- 
tion dominated; Fig. 112(1 . respectively. Thus the adiabatic 
index in this region is very close to 4/3. However, the 
fit of the physical variables by power laws as a function 
of the distance along the rotational axis produces results 
which are significantly different from those obtained for 
type-A models I3.1.1|) . and which are not in agreemen t 
with the analytic estimates of iLevinson fc Eichleil l|200n() . 
This disagreement has two reasons. On the one hand, 
the interaction between the external wind medium and 
the relativistic outflow is strongly affected by KH insta- 
bilities and not re stricted to a thin layer as assumed by 
ILevinson fc Eichlerk..200Q.) . On the other hand, the accre- 
tion torus in our models is thick (i.e., vertically extended) 
and not an infinitesimall y thin ring as in the analytic treat- 
ment bv lLevinson fc Eichler. (,2000,) . This prevents the out- 
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(a) BOl: t = 0.005 s (b) B01 : 1 - 0,010 s (c) B01 : t = 0,050 s 
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Fig. 14. Snapshots of the logarithm of the rest mass density of model BOl before the shut off of the energy release 
(the snapshot times are provided above each panel). The color scale is limited in order to enhance the details of the 
outflow region. Note that the distance scale has been adjusted in the two left panels in order to enlarge the details of 
the earlier phases of the evolution. The maximum and minimum density values (in logarithmic scale) of each model 
are given below the panels. 



flow from being sufficiently smooth to be well fitted by a 
single power law. 

3.2.1. Dependence on the total energy deposition rate 

E 

To study the dependence of the results for models of type- 
B on the assumed total energy deposition rate we have 
fixed the half-opening angle of the deposition cone to 45°. 
The models involved in this study are BOl, B04 and B05 
(Tablenj with energy deposition rates of 2 • 10^°ergs~^, 
10^^ ergs^^and lO^^ergs"^, respectively. 

Increasing the energy deposition rate in models of 
type-B leads, on the one hand, to a progressive widen- 
ing of the ultrarelativistic outflow and, on the other hand, 
to an increase of the average Lorentz factor of the fire- 
ball at all times fFig. ll5|l . From Table^we learn that the 
opening angle of the wind is smaller than the deposition 
angle ^o- Thus, as for models of type- A, the opening angle 
is not set by our choice of 0o, but is constrained by the 
presence of the thick torus. 



Within an outflow crossing time of the torus in the ax- 
ial direction (« 0.5 ms; note that the torus is more com- 
pact than in models of type- A, i.e., the crossing time is 
smaller), the fluid accelerates up to Lorentz factors of « 3, 
i.e., the Lorentz factor exceeds the inverse of the release 
angle {6^^^ = 1.27). Hence, the sideways expansion of the 
outflow is drastically reduced and the initially imprinted 
opening angle remains almost unchanged during the later 
propagation of the fireball. Indeed, if the outflow is not af- 
fected by strong KH instabilities (as in the case of model 
B04), the fluid quickly enters an ultrarelativistic regime 
where, after 10 ms, the maximum Lorentz factor can be 
larger than 30 (Table HI). 

The amount of mass in the fireball (Table ^ increases 
linearly with E. A linear fit (above the critical thresh- 
old) yields logM/ = 0.97 log - 24.2, where E and M/ 
are given in cgs units. The mass mainly accumulatess in 
the cocoon which surrounds the outflowing fireball, and 
near the leading edge of the outflow. In contrast to the 
wind solutions found for type- A models, the density varies 
by up to 10 orders of magnitude over radial scales of 
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Fig. 15. Logarithm of the Lorentz factor of the models of type-B after 0.1 s of their evolution. The color scale is limited 
in order to enhance the details of the outflow region. 

« 10* cm -10^ cm within the central cone of the fireball 3.2.2. Dependence on the half-opening angle of the 
(Fig.El)- energy deposition cone 

For a fixed, moderate energy deposition rate of 2 • 
Iq50 erg/s, we study the influence of an increasing opening 
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t = 0.100 s 




Fig. 16. Logarithm of the rest-mass density (in gem of the models of type-B after 0.1s of their evolution. The 
color scale is limited in order to enhance the details of the outflow region. 



half-angle of the deposition cone from 45° to 75° (models 
BOl, B02 and BOS; Table [2 Fig. [13). We point out that 
these models also differ in the energy deposition rate per 
unit of volume (E/Vdcp)- 



Depending on the opening angle, either a highly rela- 
tivistic outflow is produced (models BOl and B02 display 
Tmax ^ 240), or only a moderately relativistic jet is gen- 
erated (model BOS, T > 10; Fig.HSl. Therefore, the value 
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of E/Vdep determines whether a cohimated jet (for small 
values) or a wide-angle jet (for larger values) is produced. 

Type-B models display a moderately relativistic 
(r ^ 10) shell surrounding the central ultrarelativistic 
cone of the outflow, which makes it difficult to define an 
opening angle of the outflow. We define as the angle of 
the cone around the symmetry axis where F > Fop = 10. 
Using values of Fop in the interval [10 — 25] does not lead 
to substantially different values of 6w (except for model 
BOS where the maximum value of F is Fmax = 17.22). 
With this definition, and for a fixed, moderate total en- 
ergy deposition rate the opening angle of the outflow 9^ is 
sensitive to an increase of the deposition angle Oq (Table^ 
because the layer of the outflow between w 17° and « 35° 
enclosing the central relativistic or ultrarelativistic cone 
is mostly filled with matter which was pulled up and ac- 
celerated from that part of the torus overlapping with the 
deposition region. In the three models under consideration 
the deposition region extends beyond 35°. The energy per 
unit of mass released in the above mentioned layer is suffi- 
cient to speed up the fluid beyond F ^ 100, however only 
in models BOl and B02. 

After 100 ms, 6*^ ranges from 2.3° (model B03) to 30° 
(model BOl), the values decreasing with time (Tabled]). 
In models BOl and B02 this decrease is caused by the 
change from a relativistic to an ultrarelativistic flow due to 
the continuing deposition of energy. When the fluid speed 
in the radial direction approaches the speed of light the 
transverse velocity component can only be very small, i.e., 
the sideways expansion is largely suppressed. In case of 
B03, a highly coUimated jet forms very early (in less than 
20 ms), and after this wind-to-jet transition the opening 
angle of the outflow remains almost constant. Different 
from type-A models the transition is not caused by the in- 
ertial confinement due to the external medium, but by the 
growth of KH instabilities which tend to pinch the conical 
channel. These KH instabilities are of much larger ampli- 
tude for models with large because the over-pressure 
generated by the release of energy is smaller, i.e., more 
torus matter is entrained into the outflow. Indeed, the en- 
ergy in the ultrarelativistic outflow and the mass of the 
fireball AI f tends to decrease with increasing Oq. Since the 
entrained mass is not homogeneously distributed but lo- 
calized towards the lateral edges of the fireball, in these 
regions the Lorentz factor of the outflow decreases below 
10. Therefore, they are not included in our computation of 
the fireball mass (Table The process of mass entrain- 
ment is extreme in the case of model B03, where for most 
of the outflow mass F < 10 holds. 

Type-B models produce radiation dominated, wide- 
angle, shocked jets (see the profiles of the density, pres- 
sure and Lorentz factor in Fig. I12|l . Since they are ra- 
diation dominated the pressure variations effectively de- 
couple from the density fluctuations (Aloy et al. 2002). 
While the radiation pressure displays rather small vari- 
ability across and along the fireball and falls almost mono- 
tonically about 15 orders of magnitude in the radial direc- 
tion, the density variation in the same direction is more 



than 50 orders of magnitude and clearly non-monotonic 
(Fig. I12|l . The temperature of the fireball behaves similar 
as the pressure but with even less variability. The Lorentz 
factor, which also shows a non-monotonic increase with 
distance, is modulated by a number of shocks whose shape 
is either biconical (up to w 10^ cm) or radial (at larger 
distances). As a result of the non-monotonic behaviour, 
the outflow of type-B models cannot be properly fitted by 
power laws as in case of type-A models. The cause of this 
difference is easy to understand. The outflow of type-B 
models is much less dense than that of type-A models. 
Thereby, the entrainment of high density matter from the 
torus produces larger relative density perturbations than 
in type-A models. 

The differences between models of type-B and type- 
A (an d also the theoretical models of iLevinson fc EichleJ 
|2000|) are particularly large in the radially outer part 
of the fireball, while in the unshocked region the two 
types of models differ much less. This is expected because, 
on the one hand, the shocked region is neither included 
in simplified theoretical considerations nor can it be fit- 
ted by power-laws. On the other hand, the exact scal- 
ing laws of the outflow depend strongly on the details of 
the interaction between the surface of the fireball and the 
external m edium. While this interaction is assumed by 
iLevinson fc Eichlcr ( 2Qo3) to be confined to a thin layer, 
we find that it takes place in an angularly extended region 
where large-amplitude KH instabilities occur. 

As in case of type-A ni o dels ( see Sect. 13.1. 2|l the the- 
ory of ILevinson fc EichleJ l)2000() is also not suitable to 
describe the collimation of the ultrarelativistic outflow ob- 
tained in our type-B models, where the presence of the 
torus is responsible for the jet collimation, too. This fact 
is supported by the results shown in Fig. 1171 where the 
values of the half-opening angle (a) of the ultrarelativistic 
outflow are displayed calculated according to the simpli- 
fied analytic model of ILevinson fc EichleJ l|200(l) . which 
assumes that an external, baryonic wind is responsible 
for the collimation of the jet: a ~ 7r/3~^Ffj^^(L//Lu,), 
where Pn, = Vnj/c, v^, is the velocity of the wind, Ffo 
is the Lorentz factor at the base of the jet (Ffo w 1 
in all our simulations), and Lf and are the lumi- 
nosities of the fireball and of the external wind, re- 
spectively. These luminosities are evaluated as L/ = 
27ri?2c3(l - 2Rg/ Rf^'^ JH'"^'^^ d9 sin^ epVPrihV ~ 1) and 
Lu, = 27ri?2c3(l - 2i?g/i?)3/2 JJ'j^^^ de sin^ 6ipF/5,.(/iF - 1), 
respectively. In these expressions, h is the specific enthalpy 
and Pr = Vr/c with Vr being the radial velocity. For this 
quantitative comparison we identify the moderately rel- 
ativistic outflow with F < 10 around the fireball with 
the baryonic wind which produces the jet collimation in 
Levinson & Eichler's model. We further assume that the 
velocity of the wind can be approximated at each ra- 
dial distance R by the angularly-averaged, radial velocity 
of the external medium. The latter velocity is computed 
as the arithmetic mean of the radial velocities of all cells 
at a distance R in the interval between 9w{R) (the polar 
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Fig. 17. The top panel shows the luminosity of the ul- 
trarelativistic fireball (solid line) and of the surround- 
ing moderately relativistic wind (dashed line) in units of 
10^°ergs~^, as well as the angularly-averaged, radial ve- 
locity of the wind v^jc (dashed-dotted line) as a func- 
tion of the radial distance for model BOl sX t — 0.1s. 
The middle panel shows the expected fireball half-opening 
angle (in degrees) as pred icted by the analytic model of 
iLevinson h Eichleil l)2000|) . i.e., a = 7r^-ir-^(L//L^), 
where /3u, = v^jc is the velocity of the surrounding wind 
medium (w„) in units of the speed of light (c), Ffo is the 
Lorentz factor at the base of the jet (which is of order 
unity in all our simulations), and Lj and L^^ are the lu- 
minosities of the fireball and of the external wind, as de- 
fined in the text. For the calculation of a, the angularly- 
averaged, radial velocity Uu, of the external wind medium 
(see top panel) is computed at all radii R. Finally, the 
bottom panel shows |q;|° '^^ as a function of the radius 
(solid line) compared with the local half-opening angle of 
the fireball Bj as a function of the radius (dashed line). 
The value of Qj is computed as the maximum polar angle 
where, at a distance i?, the Lorentz factor of the fluid is 
larger than 10. Note that the fireball half-opening angle 
dw (given in Tables ^ and |2Jl is defined as the maximum 
of all the local half-opening angles d^. 



angle at which the Lorentz factor becomes smaller than 
10 at a radial distance i?) and 7r/2. 

Comparing the values of the fireball half-opening an- 
gle shown in Tabled with the huge (meaningless) values 



shown in the bottom panel of Fig. El we conclude that 
the non-relativistic wind is clearly not responsible for the 
fireball collimation. We point out that the value of a is 
not unique. It changes by a factor of ~ 3 depending on 
the radial distance at which it is computed. A similar vari- 
ability is also found for the local half-opening angle of the 
fireball, 0^ , calculated as the maximum polar angle where, 
at a distance i?, the Lorentz factor of the fluid is larger 
than 10 (Fig. El bottom panel). We have found that the 
function |a|° "^^ reasonably fits Qj. 

3.2.3. Dependence on the deposition volume Vjep and 
on the total energy deposition rate E for 
i?/Vdep=constant. 

In two models of type-B we have varied the half-opening 
angle of the energy deposition cone simultaneously with 
the total energy deposition rate while keeping the rate 
of energy deposition per unit of volume -E/Vdcp constant. 
The two corresponding models are B02 (with a large en- 
ergy deposition rate and deposition angle = 60°) and 
B06 (with a smaller deposition rate and — 41.4°). 

Both models develop a wide-angle ultrarelativistic jet 
because the lateral edges of the deposition region and, fur- 
ther downstream, also the lateral boundaries of the out- 
flow are causally connected, i.e., KH modes can grow mod- 
ulating the variation of the physical variables inside the 
outflow. Both models show similar global properties (prop- 
agation speed of the fireball, opening angle of the wind, 
mass of the fireball) during the time when energy is de- 
posited ( Table 0. A more detailed comparison reveals mi- 
nor differences, for example, in the maximum Lorentz fac- 
tor (larger in model B02), the propagation velocity (larger 
in B02), and the fireball opening angle after a canonical 
time of 100 ms (smaller in B02 than in BOG, although in 
both cases 6^ ~ 22°; see Tabl(^. The model with the 
smaller deposition angle (BOG) sweeps up less mass from 
the external medium, and drags less mass off the torus 
(the deposition region overlaps with the torus by w 1G°) 
than the model with the larger deposition angle and depo- 
sition rate (the energy deposition cone overlaps with the 
torus by « 35° in B02). This explains why the outflow 
in the latter model is more massive (M^°^ — 2.5 • 10^^ g) 
than in model BOG (M^o^ ^ ■ lO^^g). 

Looking at the density distribution (Fig. I1G|I we see 
that the ultrarelativistic region of model B02 is denser 
(on average) and displays much smaller density fluctua- 
tions than model BOG (Fig. El top panel). Furthermore, 
model BOG develops a more extended unshocked region 
(Fig. 112(1 , where the Lorentz factor is larger than that of 
model B02 after an evolutionary time of 100 ms (Fig. El 
bottom panel). The propagation velocity, although very 
similar, is systematically smaller in model BOG. The differ- 
ences between the two models arise from the larger angu- 
lar extent of the deposition region in model B02. Since the 
deposition region overlaps more with the torus in model 
B02 than in model BOG there is more mass pulled up from 
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the torus in the former model. This explains the larger 
density and, as a result, the smaller fluid Lorentz factor 
in the shocked region of the fireball. 

3.2.4. Dependence on the time variability of the 
energy source 

In models B07 and BOS we have varied the energy deposi- 
tion rate according to the prescriptions given in Eqs. (|10|l 
and Hll(l . The energy released in models B07 and BOS until 
time infinity is the same as in model BOl during 100 ms. 
The width of the deposition region is also equal to that 
of model BOl. With both models we try to determine the 
consequences of an intrinsic source variability as compared 
to the variability caused by the interaction with the ac- 
cretion torus and the merger halo. 

After 100 ms of continuous deposition of energy, the 
decay of the energy deposition rate of model B07 for 
t > 0.03 s has not produced any obvious signature in the 
Lorentz factor of the outflow, which in fact looks similar 
as in model BOl (Fis. I15|l. Even the variations of Lorentz 
factor in the outflow are almost indistinguishable between 
the model with constant energy release and the models 
with a variable, burst-like deposition of energy. However, 
the ultrarelativistic central core of the outflow of model 
B07 is (on average) ~ 50% denser than that of model BOl 
fFig.ll6|l. which is a consequence of a larger mass entrain- 
ment resulting from the smaller E/Vdcp near the system 
axis after ^ 35 ms relative to model BOl. 

The periodic variation of the energy deposition rate 
yields a modulated growth of the Lorentz factor of model 
BOS the periodicity following that of the energy release 
(in particular up to distances of ~ 2 • 10^ cm, where the 
grid resolution starts to become too coarse, and along 
the boundaries of the ultrarelativistic core of the fireball; 
Fig. 113). Model BOS also has a lar ger mass entrainment 
than model BOl for the reason mentioned above. It is 
slightly larger than in model B07 because of the sinusoidal 
variation of the deposition rate: E/Vdcp of model BOS is 
50% smaller than that of model B07 during the epochs of 
minimum energy release (i.e., when sin(27rt/tvar) = ^l)i 
when a larger amount of mass is entrained. Once that mass 
is inside of the fireball, it is never expelled again. 

4. Evolution after the switch-off of the central 
energy source 

In the previous sections we have discussed the evolution of 
the relativistic outflow that forms by the release of thermal 
energy around a black hole-accretion torus system. The 
simulations were carried on for a canonical time of 100 ms 
during which the source was assumed to release energy ei- 
ther at a constant rate (models BOl to B06) or with a rate 
decaying with time (models B07 and BOS) as t~^^^ (see 
? 13.2(1 . After 100 ms, the different fireballs have reached a 
distance of ~ 3 • 10^ cm. However, according to the stan- 
dard model, the GRB phenomenon can be produced only 
at radial distances of ^ 10^^ cm from the energy source. 



i.e., the outflow must still propagate and persist over an- 
other 4 orders of magnitude in radius before an observable 
signal emerges. In order to investigate longer timescales 
than the epoch of energy release, we have followed the 
evolution of the type-B models for another 400 ms after 
the central energy supply had either been shut off (BOl 
to B06) or during which it continued at a decreasing rate 
(models B07 and BOS). 

In the following sections we will concentrate on the 
evolution of the prototype model BOl I4.1|l . and after- 
wards we will outline the particularities of the other type- 
B models f? 0|l . 

4.1. Model BOl 

After the epoch of 100 ms of a constant rate of energy re- 
lease has ended (except for the first 10 ms during which 
the energy deposition rate is linearly increased) , the high- 
Lorentz factor outflow possesses a conical shape which 
decouples from the energy deposition region, and moves 
away from the black hole (Fig. 118(1 . The shape is mainly 
a consequence of the continuing radial acceleration of the 
matter in the outflow. The angular structure of the adi- 
abatically expanding fireball is non-uniform because of a 
larger baryon loading at larger polar angles. The Lorentz 
factor reaches values between 20 and 50 (approximately, 
between « 15° and 30°) where the fireball is loaded with 
baryons extracted from the outermost layers of the torus. 
These interaction layers have still mildly relativistic or 
subrelativistic velocities at the time the energy release 
is switched off. In contrast, the ultrarelativistic core has 
swept up only matter from the dilute halo. One, there- 
fore, finds variations of the Lorentz factor fFig. llS(l which 
depend on the lateral angle (in addition to the radial vari- 
ations of r). The maximum Lorentz factor (Fmax — S50) 
is attained close to the lateral edges of the fireball core (at 
polar angles around S°) and in its rear part, while in the 
region around the axis {9 < 3°) the fluid Lorentz factor is 
about 20% smaller than its maximum value fFig. 121(1 . For 
small angles, the Lorentz factor exhibits a fast rise from 
the backward end of the flreball to its center where values 
beyond 800 are found fFig l20(l . From this local maximum 
the Lorentz factor decreases slightly towards the terminal 
shock of the wind (see the violet shade in the right panel 
of Fig.[T5)l. For larger angles {9 > 15°) the radial profile of 
the Lorentz factor is similar although the values here are 
much smaller (F ~ 30). This radial profile of the Lorentz 
factor is the result of the beginning of a coasting phase, 
that starts soon after the energy release ceases. 

The rest- mass density (Fig. 119(1 in the fireball does not 
preserve the irregular structure imprinted by the KH in- 
stabilities that had developed near the lateral edges dur- 
ing the earlier evolution. After the shutdown of the en- 
ergy deposition, two voids (density minima) show up in 
the central part of the fireball, one close to its rear part 
and another one close to the front of the fireball (solid line 
in FigEOl- The two voids are radially separated by a 20 
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(a) BOl: t = 0.100 s (b) B01 : t - 0,£50 s (c) B01 : t = 0.500 s 
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Fig. 18. Evolution of logarithm of the Lorentz factor of model BOl after the shut down of the central energy supply. 
Time is measured from the moment when the energy deposition was started. 



times denser "wall" , corresponding to a reverse shock that 
sweeps backwards into the fireball. 

The pressure in the fireball decreases with increasing 
radial distance and increasing polar angle. Across the fire- 
ball, in the polar direction, the pressure is lowest close to 
the side edges of the fireball core and is highest around 
the symmetry axis (where the density is about a factor 
of 2 larger than at the sides; Fig. I21|l . The pressure vari- 
ation in the radial direction can be as large as three or- 
ders of magnitude (Fig|2DJ). However, due to the fact that 
the fireball is radiation dominated, the variation of pres- 
sure (or temperature) is much smaller than that of the 
rest-mass density. The structure of the outflow can be 
understood by the combination of three main dynamical 
effects. First, there is a trend to adopt a radially strati- 
fied wind structure due to the expansion of the outfiow as 
it moves radially outwards. This explains the decrease of 
the pressure, temperature and density (and the increase of 
the Lorentz factor) in the radial direction. Second, there 
are local variations of the physical variables caused by 
the KH instabilities during the early interaction of the 
fireball with the external medium and, particularly, with 
the torus. The largest of these instabilities survive beyond 
the shutdown of the energy deposition. Most of the small 



scale variations are (numerically) suppressed because of 
the logarithmic coarsening of the grid in the radial direc- 
tion. However, KH modes with wavelengths smaller than 
the width of the thick, expanding layer surrounding the 
central ultrarelati vistic cone are also physically damped 
l(Birkinsh a-^|l|^gl [). Third, there is a genuinely rela tivis- 
tic effect ((Pons et al."2000: lR,ezolla fc Zanott i' '2002*) that 
leads to the generation of a rarefaction wave and a shock 
(instead of two shocks) for collisions between fluids in 
case of a small impact angle (an example of this phe- 
no menon in the con text of parsec scale jets can be found 
m lAlov et ablEnnl . The effect arises in the presence of 
a non-vanishing tangential velocity component at a rel- 
ativistic shock. The generalized jump conditions require 
that the quantity hTvt is equal (where h and vt are the 
specific enthalpy and the velocity component tangential to 
the shock front, respectively) on both sides of the shock. 
The rarefaction propagates from the collision region into 
the medium that has smaller density and higher specific 
enthalpy. The shock moves in the opposite direction into 
the fiuid that has larger density and smaller specific en- 
thalpy. 

This described situation applies to our models when 
the fireball expands sideways (i.e., in the lateral direc- 
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(a) B01M = 0.100 3 (b) B01: 1 ~ 0,250 s (c) B01 : t = 0,500 s 
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Fig. 19. Evolution of the logarithm of the rest-mass density (in gem of model BOl after the shutdown of the central 
energy supply. Time is measured from the moment when the energy deposition was started. 



tion) and interacts with the external medium and the 
torus. Tangential velocity components are present due to 
the radial fluid motion in the fireball, the moderately rel- 
ativistic rotation of the torus, and the expansion of the 
nonrelativistic ejecta that are blown off the torus surface 
due to the energy deposition. The interaction leads, on 
the one hand, to the formation of a rarefaction wave that 
starts from the lateral edges of the fireball, propagates 
very slowly towards the axis, and leaves behind a region 
with negative angular (6*—) velocities (i.e., velocities di- 
rected to the axis). On the other hand, the interaction 
generates a shocked layer that encloses the fireball's side 
edges. Because of the rarefaction wave the pressure near 
the edges of the fireball core (which is a region of low 
density and high specific enthalpy) is smaller than near 
the axis (in Fig. 1211 the rarefaction extends from « 7° to 
« 14°). Analytic modeling of the collision of the fireball 
with the external medium in two dimensions (assuming 
planar symmetry) shows that a rarefaction wave forms in- 
stead of a shock when the impact angle is smaller than 
- 35°. 

We point out that the existence of such a rarefac- 
tion wave w as included neither in analytic models — 
iLevinson fc Eichler. (,200fl>) assumed that a thin shocked 



layer that results from the propagation of two oblique 
shocks — nor in phenomenological fireball-jet models as- 



omogeneous (e.g., 


Piran 


Il999ll200(* 


iKumar & Granot 


I2OO3 


; ISalmonsonl 



I2nn.1l) . 

From the above discussion it is clear that the ultra- 
relativistic outflow is laterally inhomogeneous and not 
bounded by sharp edges (as assumed in the so-called ho- 
mogeneous or top-hat jet models). The outflow profile can 
also not be fitted by a monotonic function (a power-law 
or a Gaussian function as assumed in the structured jet 
models). Let us now consider the one-dimensional angu- 
lar outflow structure by calculating the radial average of 
the quantities displayed in Fig[53 We apply the formula 
A = (AVr)-' jl^--_r^Adr, with = (i?L.-^LJ/3, 
Rmax and Rmin being the radii of the front and of the 
rear end of the relativistic outflow, respectively. We find 
that the evolution of the angular pro files is rather dif- 
ferent from the behaviour discussed bv lKumar fc GranoiJ 
l|2nn3() . The angular wedge with high Lorentz factors be- 
comes narrower with time and the angular location of the 
peak Lorentz factor moves to smaller angles. At the same 
time the peak value increases and becomes larger than the 
value on the symmetry axis (Fig. [23) . A good fit for the 
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Fig. 20. Logarithm of the rest-mass density (top panel), 
of the pressure (middle panel) and of the fluid Lorentz 
factor (bottom panel) vs radius along the symmetry axis 
for models BOl, B07 and BOS (see legend in the plot) after 
500 ms. 

radially averaged Lorentz factor (f ) as a function of the 
polar angle in the interval [0 
sion 

r = 1 



, 90°] is given by the expres- 



flo + aix + a2X + a^x 



(12) 



04 + a^x + exp (aex) 

The values of the seven coefficients Ui and the of 
the fit are listed in Table |3 for the type-B models at 
a time of 0.4 s after the shutdown of the central energy 
source. Attempts to fit simpler functions in a narrower in- 
terval (namely, [0°,20°]) were unsuccessful. In particular, 
a Gaussian profile of the form ae -1-1, or a quadratic 
ansatz of the form a/{b + cx'^) had to be rejected because 
the resulting values arc 176 and 860, respectively. 

Figure 1221 shows that during the last 0.4 s the increase 
of the averaged Lorentz factor per unit of time in a 10° 
cone around the axis is much smaller than during the first 
0.1 s. However, we also sec that the averaged Lorentz fac- 
tor has not yet converged to an asymptotic value after 
0.5 s, i.e., acceleration continues. Table|21gives a rough es- 
timate of the maximum bulk Lorentz factor (defined in the 
table caption) that can be asymptotically reached by the 
fireball when its whole internal energy is converted into 
kinetic energy (provided there is no further acceleration 
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Fig. 21. Lorentz factor (dashed line), pressure (dotted 
line) and rest-mass density (solid line) in polar direction 
at a radial distance of 1.4 ■ 10^" cm after 500 ms of evolu- 
tion of model BOl. The quantities are normalized to their 
respective maxima in the displayed angular interval. Note 
the rarefaction visible between the center and the edges 
of the fireball (at « 8°). 

or deceleration due to the presence of a pressure gradi- 
ent in the ambient medium of swept-up matter, and that 
the medium at larger distances from the fireball origin is 
uniform and very diluted). The values of Too are still ^ 2 
times larger than the corresponding values of Fmax except 
for model BOS (Table This is consistent with our hy- 
pothesis according to which the value of Poo sets an upper 
limit to the maximum bulk Lorentz factor of the outfiow. 
Furthermore, the fact that Fmax < Foo suggests that there 
is still room for additional acceleration of the fireball fluid. 

Profiles of the radially averaged energy density (e) can 
be fitted as a function of the polar angle, in the interval 
[0° , 90°] , by polynomials of the form 



4 

i=0 



(13) 



The values of the five fit coefficients hi 0.4 s after the 
shut down of the energy deposition are listed in Table 01 
Simple second-order polynomials provide good fits to e in 
the interval [0°,20°]. 

Radial averages of the angular component of the four- 
velocity (Tve) display profiles with negative values in some 
regions of the body of the fireball. This is a consequence 
of the above mentioned rarefaction wave that results from 
the interaction between the external medium and the fire- 
ball. These negative velocities tend to collimate the out- 
flow. Hence, as Tve is a measure of the lateral expan- 
sion velocity in the comoving frame, it turns out that 
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Table 3. Fit coefficients of the lateral profiles fEa. I12|l of the radially averaged Lorentz factor 0.4s after of the cease 
of the energy release. Note that model BOS cannot be matched with the functional form given in Eq. H12|l and. 
provides a measure of the quality of the fit. 
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Table 4. Fit coefficients of the lateral profiles of the ra- 
dially averaged energy density fEa. I13|l at i = 0.5 s. 
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140 


-82 


0.003 
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19 
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94 
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the part of the fireball which has average Lorentz factors 
10 ^ r ;S 100 is contracting while matter with relatively 
large Lorentz factors ( }t 100) expands with velocities that 
are considerably su bsonic (Fig. | 22|). In this p oint we qual- 
itatively agree with lKumar fc GranotI l)2003l) : The lateral 
expansion is slower than that assumed by the simplest 
homogeneous models (which consider expansion velocities 
close to the sound speed or even larger). We also point 
out that the absolute values of Tvg tend to decrease with 
time. 

An upper bound for the duration of a GRB event might 
be set by the time interval between the moments when 
the forward and the rear ends of the expanding fireball 
reach transparency at large radii. This time difference is 
linked to the difference in speed between the forward and 
rear ends of the fireball in the radial direction. It turns 
out that the velocity of propagation of the forward edge 
corresponds to that of the terminal shock of the wind or 
the head of the jet in our models. This velocity is deter- 
mined by the amount of mass swept up by the bow shock 
of the jet. Since the halo density in models of type-B is 
low, the propagation speed of the fireball is ultrarelativis- 
tic (Table 0. On the other hand, the velocity of the rear 
part of the wind is set by the speed of the contact dis- 
continuity that separates the fireball from the external 
medium after the energy release has stopped. Behind this 
wave, matter fills the wake of the fireball, i.e., the ultrarel- 
ativistic wind structure is lost. The contact discontinuity 
propagates at a speed very close to the speed of light be- 
cause the relativistic wind is very rarefied. In addition, 
the speed of this wave is not constant, because the plasma 



in the fireball experiences ongoing acceleration during the 
transformation of its internal energy into kinetic energy. 
In order estimate this velocity we have approximated the 
hydrodynamic states at the rear, radial edge of the fire- 
ball by two constant states corresponding to the fireball 
(right) and to the external medium (left). The two states 
are (in cgs units): (pL,PL,rL) « (10i'*,4 ■ 10^^,2) and 
(pr, Pr,^r) ~ (6 • 10^^, 10^'', Fr). The generic solution 
of the Riemann problem (see above) consists of a con- 
tact discontinuity flanked by two rarefaction waves (one 
advancing into the fireball and the other one extending 
into the external medium) where the fluid progressively 
accelerates from Fl to Fr. We have computed the ex- 
act solution of the Riemann problem set by this series of 
states as a function of Fr with 80 < Fr < 220, which are 
typical values of the Lorentz factor near the rear edge of 
the fireball as it accelerates after the energy deposition is 
stopped. The Lorentz factor associated with the contact 
discontinuity (Fcd) is a monotonic function of Fr that 
can be fitted by the linear relation 

Fcd « 6.42- 10"2Fr + 12.60, (14) 

i.e., Fcd < Fr for ultrarelativistic fireballs. 

Fr is, in turn, smaller than the Lorentz factor asso- 
ciated with the leading radial edge of the fireball during 
the acceleration phase of the outflow, which leads to a 
progressive stretching of the fireball in the radial direc- 
tion. A possible important consequence of this stretching 
is that the duration of the GRB event might be consid- 
erably longer than the time interval of the activity of the 
central engine. If the fireball stretches radially on its way 
to the optically thin regime, the on-time of the energy 
source would not necessary be a suitable measure of the 
GRB duration. In order to derive an order of magnitude 
estimate of the stretching effect it is necessary to extrap- 
olate from our computed evolution to much later times. 
For this purpose we have calculated from our hydrody- 
namic models the speed of the fireball tail and compared 
it with the analytic estimate of Eq. p4(l . It turns out that 
the analytic result agrees reasonably well with the numeri- 
cal computation. The contact discontinuity, which roughly 
corresponds to the tail of the fireball, moves with an av- 
erage velocity of vcd — O.OOlg'™^. Hence, if the burst is 
produced after the fireball has traveled a distance R, an 
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Fig. 22. Radially averaged values of the Lorentz factor Fig. 23. Radially averaged values of the Lorentz factor 
(top panel), angular component of the four- velocity (in- (top panel), angular component of the four- velocity (in- 
termediate panel) and energy density (bottom panel) as a termediate panel) and energy density (bottom panel) as 
function of the polar angle for model BOl. The different functions of the polar angle after 0.5s of evolution. The 
lines correspond to different moments in time spaced by different lines correspond to different models from BOl to 
100ms (the times are listed in the top panel). B06. 



upper bound to its duration will approximately be given 

by 

AU< (^R-'^(R-U,c)^c-\ (15) 

where Vp is the propagation velocity of the fireball in the 
radial direction and tga = 0.1 s is the time of source activ- 
ity. 



Taking a typical distance R = 10^^ cm, Eq. (|15|) yields 
a maximum burst duration of Atf, ~ A.S'^^^^q^ s, which is 
considerably longer than the time interval of source activ- 
ity {tsa = 0.1 s). It is emphasized that Eq. (|15|l and the 
arguments given in this context do not allow one to make 
quantitatively meaningful predictions for observable GRB 
durations. The estimate of Eq. (|15|l assumes that no de- 
celeration of the fireball head occurs by the accumulation 
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Fig. 24. Isotropic equivalent total (left) and kinetic (right) 
energies of matter with T > 100 for all B-type models 
at different evolutionary times. Model BOS is omitted in 
this figure because it does not develop an outflow with 
r > 100. Note that the central plateau of the isotropic 
equivalent total energy increases substantially between 
0.1s and 0.2 s when the source is already off. This is due 
to the continuous acceleration of matter to Lorentz fac- 
tors above 100 as a result of the conversion of internal to 
kinetic energy in the fireball (even after the shut-off of the 
energy release). 

of matter that is swept up from the ambient medium un- 
til the fireball reaches the distance R. A detailed model of 
the long-term evolution and structure of the fireball also 
requires a realistic, time-dependent model for the torus 
accretion and the associated energy release (see the sub- 
sequent discussion of models B07 and BOS). Moreover, the 
velocities of the head and the tail of the fireball cannot be 
extracted to very high accuracy from our rclativistic hy- 
drodynamics results (mainly because of the limited resolu- 



tion at large radial distances), but are nevertheless used to 
extrapolate the simulations by several orders of magnitude 
in radius. The velocities also do not need to be decisive for 
the actual duration of an observable GRB, in particular, 
because the tail velocity was deduced from the motion of 
the rear contact discontinuity of the fireball, which might 
not be the backward boundary of highly relativistic mat- 
ter with r > 100 that forms the 7-ray emitting part of 
the GRB fireball. Equation H15|l is mainly meant to mani- 
fest the fact that the GRB duration, in particular for short 
bursts, is not necessary given by the on-time of the energy 
source, but can be significantly longer due to the discussed 
fireball stretching. Indeed, the stretching may account for 
the major part of the GRB duration. The numerical es- 
timates obtained from Eq. H15|l at best are rough upper 
limits to this effect. In contrast, in the coUapsar model for 
long bursts, for example, the accretion time and thus the 
period of source activity clearly dominates the observed 
GRB duration. 

4.2. Models B02 to B06 

This series of models has in common with BOl that the 
source activity time is limited to 0.1s. In this section we 
shall discuss which are the similarities and differences of 
their evolution after the end of the energy release rela- 
tive to model BOl. In Fig. 1531 the fluid Lorentz factor of 
all type-B models is displayed at the end of our computed 
evolution (0.5 s). The high Lorentz factors suggest that all 
models, except BOS, may be suitable for producing cosmic 
GRBs (at least they fulfill the canonical requirement of 
F }t 100). All models have outfiow opening angles which 
are smaller than the opening angle of the deposition region 
(see also Table confirming once again that the opening 
angle of the outflow is set by its interaction with the side 
walls of the torus. The diversity of opening angles (from 
3° to 25°) is caused by the differences in the energy that 
is released around the BH per unit of time and unit of 
volume, which also depends on the angular width of the 
deposition region. The trend is that models with a larger 
E/Vdcp develop larger outflow opening angles. Model BOS 
is special only by the fact that due to the large opening 
angle of the deposition cone and a modest total rate of en- 
ergy release, most of the energy is deposited in the densest 
part of the torus and hence can only trigger a moderately 
relativistic outflow. In case of models B02 and B04 the 
final opening angle of the ultrarelativistic outflow, might 
be smaller than at 0.5 seconds because both models show 
substantially negative values of Tvg (« —0.5c; Fig. I2S|) . 

Stretching of the fireballs in the radial direction should 
occur in all models discussed in this section, and is indeed 
found up to 0.5 s. The radial extension increases similarly 
in all models. In agreement with our numerical results and, 
according to our estimate of Eq. H14I) , the small differences 
result from somewhat different velocities of the tail of the 
fireballs (the forward radial front of the fireballs moves in 
all the cases except BOS at Vp ~ c). 
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Fig. 25. Logarithm of the Lorentz factor of the type-B models after 0.5 s of evohition time. The color scale is limited 
in order to enhance the details of the outflow region. 



The maximum Lorentz factor found within the fireball Lorentz factor. Figure El shows that soon after the end 

a.t t = 0.5 s varies substantially from model to model and of energy release the total (internal plus kinetic) isotropic 

lies in the range from ~ 16 to ~ 1000 (Table|2l, following equivalent energy of the fireball core reaches a maximum 

the trend that a larger E/Vdcp leads to a larger maximum and saturates, while there is an ongoing transformation 
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t = O.BDO s 

Fig. 26. Logarithm of the rest-mass density (in gcm^'^) of the B-models after 0.5 s of evolution time. The color scale 
is limited in order to enhance the details of the outflow region. 

of internal into kinetic energy (which is constantly ris- A comparison of the estimated terminal Lorentz fac- 
ing). During this conversion of energies the Lorentz factors tors Too of the different models (Table |2J) reveals that 
within the fireball grow. They have not saturated until the models BOl, B04, B05 and B06 are somewhat more ef- 
end of the simulations. ficient than models B02 and BOS, the main reason being 
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the larger opening angle of the energy deposition cone in 
the latter two cases in combination with a smaller energy 
deposition rate per volume. This causes an increase of the 
baryon loading of the fireball relative to the energy stored 
in the pair-photon plasma and is a crucial difference be- 
tween models BOS and, e.g., B04. The latter develops an 
energetic ultrarelativistic jet although the energy deposi- 
tion per solid angle is smaller by a factor of ^ 10 than in 
BOS. In model BOS the estimate of the asymptotic Lorentz 
factor is not meaningful because the model develops a col- 
limated jet with only a very narrow and light relativistic 
core (where F ^ 10; Fig. l2S|) . As the estimate Too only con- 
siders the mass of the outflow with F > 10, the amount 
of mass that fulfills this criterion in model BOS is dynami- 
cally negligible compared to the mass of the whole outflow, 
which would have to be used for calculating the asymp- 
totic value of F in this case. 

All models (except BOS) display a similar structure 
of the radially averaged energy density (Fig. I2S(I . Close 
to the symmetry axis and up to a few degrees off axis 
it is nearly constant. At larger angles 9 a hump appears 
where the energy density is about a factor of three higher 
than near the axis. From this local maximum the energy 
density declines gently towards the lateral edges of the 
fireball where it starts rising again due to the interaction 
with the external medium. The energy hump within the 
fireball coincides with a region of large compression of the 
fluid where the comoving lateral velocity reaches a local 
minimum that, in some models (e.g., in B02, B04 and 
BOS), is an absolute minimum; (Fig. ESI middle panel). 

The isotropic equivalent total energy as a function of 
polar angle 6 for matter with Lorentz factors larger than 
50 (Fig. I24|l consists of two parts: a central core of 3° 
(model B02) to 9° (model BOl) where it is nearly uniform 
and maximal, surrounded by an extended layer where it 
drops roughly like 9~^. Considering that the equivalent 
isotropic total energy is proportional to the energy per 
unit of solid angle, our energetic profiles are in very good 
agreement with those assumed in the st ructured jet model s 
for the afterglows of long GRBs (e.g., ISalmons'oi]l200Sl) . 
However, we point out that our results are supposed to 
apply for the GRB event itself and for short bursts. The 
generation of the energy profile of the burst seems nei- 
ther to depend on the width of the deposition region (it 
is similar for a wide range of Oq values) nor on the abso- 
lute value of the energy deposition rate (or its time de- 
pendence). Nevertheless, we are not in a condition yet to 
claim that this is a universal profile, because it might be 
linked to the particular spatial dependence of the assumed 
energy deposition law or to the properties of the BH-torus 
system. Future simulations will have to investigate these 
aspects. 



4.2.1. Models B07 and BOS 

In these models, the energy is released with a time- 
variable deposition rate fEa [T?Hl and continues for longer 
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Fig. 27. Radially averaged values of the Lorentz factor 
(top panel), angular component of the four- velocity (in- 
termediate panel) and energy density (bottom panel) as 
functions of the polar angle after 0.5 s of evolution. The 
different lines correspond to models BOl, B07 and BOS, 
respectively. 



that 100 ms. Note that due to the assumed power-law 
decay the energy deposition rate is still at a level of 
3.5 • 10^^ergs~^after 500ms. This is the reason why the 
high-Lorentz factor outflow that has formed an almost 
conical fireball docs not detach from the energy deposition 
region and, thus, the axial funnel of low density and high 
Lorentz factor stays open until the end of the computed 
evolution fFigs. 1^ I29I2F|| . This axial spine decreases 
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in width with time and, if we continued to foUow the evo- 
lution of B07 or BOS for a longer period, the energy supply 
at some later time would not be sufficient any longer to 
keep the funnel open. 

The structure of the fireball is non-uniform 
(Figs. EHl and its morphology does not differ 

substantially from that of model BOl (except for the long 
trailing tail which results from the continuous energy 
deposition). The plateau of the maximal, radially aver- 
aged Lorentz factor is more narrow (w 2° in both model 
B07 and BOS) than in model BOl (Fig. because of 
the progressive contraction of the fireball fluid (which is 
larger in model B07 than in BOS; Fig. 1271) . Furthermore, 
the radially averaged maximum Lorentz factor after 0.5 s 
is larger for these two models (with values above 375) 
than for BOl. However, the absolute maximum Lorentz 
factor (Fniax) at the end of the simulations corresponds 
to model BOl (Table EJ- The reason for these differences 
is that BOl is more efficient than B07 or BOS in using the 
released energy to accelerate matter to Lorentz factors 
above 100 (see the values of Eryioo/Ed in TableEJ. 

The estimates for the asymptotic Lorentz factor Too 
for models B07 and BOB in Tabic |2 may not be meaning- 
ful because energy is released even beyond the simulated 
epoch of evolution. 

The time variability of the energy deposition rate of 
model BOS produces a wrinkled boundary of the fireball 
with a modulation period of 10 ms as the energy release 
(Fig.|5Sl left panel). A similar wrinkled surface occurs in 
model B05 (Fig. I16|l . but different from model BOS the 
modulation of the density in the surface is not exactly pe- 
riodic, and the values of the density in the fireball core up 
to ~ 2 • 10^ cm are much higher 10~^gcm~'^ very near 
the axis; Fig.[lBl). After 0.25 s models B07 and BOS display 
a fairly similar structure. This is, to a large extent, due to 
the grid coarsening in radial direction. Beyond a distance 
of « 10^" cm the grid spacing exceeds 10 light-milliseconds 
and, therefore, most of the variability on scales less than 
« 15 ms is damped or erased. For this reason, we can- 
not draw firm conclusions on the short-time variability 
from our models, and in particular we cannot answer the 
question how it depends on the intrinsic variability of the 
source or the interaction with the external medium. 

The radially averaged energy density of models BOl, 
B07 and BOS varies somewhat differently with the polar 
angle (Fig. 127(1 . The central, most energetic core of the 
fireball of models B07 and BOS is more narrow. Model 
BOS does not possess a clear plateau around the symmetry 
axis and also has, on average, a smaller energy content 
compared to the other two models. The kinetic energy in 
the central core of the ultrarelativistic outfiow is smaller 
than in the surrounding layer in all B-models at t = 0.5 s 
(Fig. because the higher Lorentz factors in the cores 
are compensated by larger densities outside. The efficiency 
for converting thermal to kinetic energy until ~ 0.5 s is 
^ 25% smaller in case of model B07 than in model BOl 
or BOS (Table[2Il. However, in all three cases and all other 
B-models the kinetic energy in the fireball relative to the 



total energy of matter with F > 100 is still at most a few 
per cent at t — 0.5s (Table I^J, i.e., there is ample room 
for further fiuid acceleration. Figure 1301 shows that only 
model BOS is an exception, as it has essentially reached its 
terminal state. Figure 1201 also shows that the contribution 
to the kinetic energy of the parts of the outflow having 
bulk Lorentz factors larger than 50 is small while most 
of the kinetic energy is transported by fluid with Lorentz 
factors below 10. It is also evident that the fraction of 
kinetic to internal energy at 0.5 s is closer to unity for 
matter with low Lorentz factors farther off axis, i.e., this 
material will not reach much higher asymptotic velocities 
during its adiabatic expansion. The ultrarelativistic core, 
however, can be expected to continue accelerating, when 
its internal energy is converted to kinetic energy. 

The evolution of the ratio of kinetic energy to total (in- 
ternal plus kinetic) energy in the outfiow of the prototype 
model BOl is shown in Fig. |^ In this figure the kinetic 
and total energies of the relativistic outfiow are computed 
as the integral values for all cells with F > Fmin for differ- 
ent lower limits Fmin of the Lorentz factor. One can notice 
a slight, non-monotonic increase of the fraction of kinetic 
energy with time. Several conclusions can be drawn from 
Fig. EH First, if the rate of increase can be extrapolated 
until the time when transparency sets in, it must be ex- 
pected that most of the internal energy will be converted 
to kinetic energy. Second, after 0.5 s of evolution about 
half of the kinetic energy of the outflow is in fiuid with 
r > 10. However, it is impossible to reliably predict the 
final distribution of kinetic energy versus Lorentz factor 
(after, e.g., ~ 1000s), because the trends shown for fiuid 
elements with different values of Tnun are slightly different 
and, in some cases it is hardly possible to deduce an ac- 
curate number for the net increase of the ratio of kinetic 
to total energy (e.g., in case of Fmin > 100; Fig. rTT|l . 

5. Summary and conclusions 

We have performed general relativistic, axisymmetric hy- 
drodynamic simulations for studying the creation and 
evolution of relativistic outflow in response to the de- 
position of thermal energy above the poles of a black 
hole-accretion torus system. The latter is considered to 
be a necessary condition for obtaining baryon-poor fire- 
balls, in particular if the energy is provided by neutrino- 
antincutrino annihilation. Above the poles of the black 
hole a favorable environment is encountered for the cre- 
ation of an ultrarelativistic e^-pair-photon plasma by the 
mentioned process, different from the situation in case 
of neutron stars, where neutrino energy deposition in 
the surfa ce-near lay ers produces a nonrclativistic bary onic 
wind tPuncan et al.-19 S6:■^Yoosley 1993;.Qian fc WooslevI 
hflflftlThomnson et a].ll20(m 

Be ing g uided by numerical models of lRuffert fc Jankal 
1)1999(1 and lJankaet al., (1999) , the BH-torus systems for 
our studies were set up to represent the remnants of 
NS-I-NS or NS-I-BH mergers. We chose Schwarzschild BHs 
of 2.44 Mq or 3 Mq, being girded by a thick torus of 
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(c) B08n:t = 0.500 s 
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Fig. 28. Evolution of the logarithm of the rest-mass density (gem ^) of model BOS after 100ms. Time is measured 
from the moment when the energy deposition was started. 



0.13 Mq or 0.17 Mq, respectively. The accretion torus 
was constructed either by relaxation of an initial toroidal 
cloud of matter or by creating a near-equilibrium con - 
figuration with the technique of iFont fc Daign3 l|2002|) . 
In the first case the BH-torus system is surrounded by 
a fairly dense (~ lOOgcm"^) gas halo with a radius of 
about 10^ cm (containing a mass of roughly 10~'*Mo), 
e.g., considered to be the consequence of a neutrino- 
driven wind from a transient, metastable NS which un- 
derwent collapse to a BH-torus system only after some 
time delay. In the second case the density of the sur- 
rounding gas is much lower and decreases radially ap- 
proximately as r~^'^. Correspondingly, we discriminate 
between high-density type-A models and low-density type- 
B models. Both are considered as possible results of a vi- 
olent, preceding merger history. Thermal energy was de- 
posited within a cone with chosen opening angle above the 
poles of the black hole, starting at the inner grid boundary 
at 4G'Mbh/c^ or 2GMbh/c^ for type-A or type-B mod- 
els, respectively. Cases with constant energy deposition 
rate were studied or with a burst-like initial phase, fol- 
lowed by a long-time, gradual decay. The energy depo- 
sition rate was assumed to depend only on the vertical 
height z measured along the rotation axis and to drop 



like . This prescription is an approximate representa- 
tion of deposition maps obtained by detailed calculations 
of ncutrino-antineutri no annihilation in the vicinity of 
the BH-torus system dRuffert &: Jankalll999t Ijanka et alJ 
Il999) . Disregarding the modest lateral bending of the 
computed surfaces of constant vv-eTieigy deposition rate 
turned out not to be essential for the discussion in this 
work, but the deviation of these surfaces from spheres had 
important consequences. 

We calculated sequences of type-A and type-B mod- 
els, systematically varying the total rate of energy deposi- 
tion [E) while keeping the lateral width of the deposition 
region unchanged. Alternatively, we changed the opening 
angle of the energy-deposition cone for a fixed integral rate 
or for a fixed energy deposition rate per unit of volume 
(i?/Vdop), respectively. In summary, our studies revealed 
the following dynamical behavior of the evolving outflow: 

— Relativistic outflow develops in form of narrow, knotty 
jets or wide, ultrarelativistic winds or jets. The gener- 
ation of relativistic outflow requires a high energy de- 
position rate or low environmental density. In case of 
the high-density type-A models the energy deposition 
rate has to exceed a certain threshold value (around 
lO^^ergs"^) in order to produce any outflow. Slightly 
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Fig. 29. Evolution of the logarithm of the Lorentz factor of model B07 after 100 ms. Time is measured from the 
moment when the energy deposition was started. 



above this threshold an inhomogeneous jet is formed 
rather than an ultrarelativistic wind. 
In case of type-A models the collimation of the outflow 
is mainly determined by the interaction with the high- 
density ambient medium. The opening angle of the 
polar outflow increases with E/Vdcp and varies from 
a few degrees for the low-density (baryon-poor) jets 
that develop just above the ejection threshold to more 
than 25° for the smooth, low-density, wide-angle winds 
in case of the highest considered E/Vdcp- The interac- 
tion between outflow and ambient medium resembles 
the hydrodynamic collima tion mechanism discussed by 
iLevinson fc Eichlerl ^0(f), and the unshocked part of 
the baryon-poor outflow can be well fitted by power 
laws with indices close to the ones predicted by the 
analytic model of these authors. Because of the large 
mass that is swept up, the propagation speed of the 
shocked, forward part of the outflow becomes only 
mildly relativistic in all models, not allowing for fa- 
vorable conditions to produce GRB emission. Instead, 
a low-luminosity (at most ^ lO'^'^ergs"^) soft UV- 
flash at a temperature of ~ 5 • 10^ K with a duration of 
~ 1000 s can be expected as the observable signature of 
merger events surrounded by high-density halos. Due 



to their small luminosity, only galactic events might be 
detected. 

In case of type-B models all employed energy deposi- 
tion rates (above ~ 10**® erg s"""^) led to non-uniform, 
wide-angle, radially and laterally structured ultrarela- 
tivistic jets provided the rate of energy deposition per 
unit of volume was sufficiently high. The propagation 
velocities are found to be very close to the speed of 
light, corresponding to Lorentz factors of several hun- 
dred, also in the forward, shocked part of the flow. 
This suggests favorable conditions for the production 
of GRBs . The simple power - law f its of the analytic 
model of iLevinson fc Eichleil ( 200(1) do not work well. 
The latter authors assumed that a baryon-poor jet is 
confined by the ram pressure of a baryonic (nonrel- 
ativistic) wind which originates from a toroidal ring 
of negligible thickness. They further assumed that the 
baryon-rich wind collides with the relativistic jet and 
gets deflected into a very thin, shocked layer that en- 
velopes the relativistic flow. Our hydrodynamic simu- 
lations do not support this picture. We found that the 
collimation of the baryon-poor flow in type-B models 
is caused by the interaction of the accelerating fluid 
with the dense accretion torus, whose vertical exten- 



Aloy, Janka and Miiller: Production of short GRBs 



37 




% m m M % so ^ ^ 



Fig. 30. Isotropic equivalent total energy (left panels) and isotropic equivalent kinetic energy (right panels) correspond- 
ing to outflows with Lorentz factors above certain thresholds at 500 ms of computed evolution. Rcsuls for different 
thresholds are plotted with different line styles. Note that in the panels corresponding to model BOS there are no lines 
for Lorentz factor thresholds above 50, because there is no outflow with such a large Lorentz factor. 



sion cannot be ignored. This interaction produces a 
thick shell of gas which surrounds the ultrarelativis- 
tic core and which shows large inhomogeneities due to 
Kelvin- Helmholtz instabilities. The half-opening angle 



of the ultrarelativistic flow is set within the first mil- 
lisecond of its expansion on a distance that corresponds 
to the vertical thickness of the accretion torus. It is 
determined by the curvature of the torus walls around 
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the poles of the BH and was found to be around 20° 
in the computed type-B models. 
— Type-A and type-B models also showed a different be- 
havior after the deposition of energy was switched off 
or had decayed. In case of a high-density halo around 
the BH-torus system (type-A) the conical, relativistic 
wind structure is destroyed within a few source activity 
timescales (tsa), because the tail of the outflow moves 
into a cleared funnel and is much faster. Therefore, it 
catches up with the leading edge of the flow, which gets 
decelerated by swept-up matter. In contrast, in type-B 
models with their steeply decreasing halo density, the 
front part of the outflow is steadily accelerated while 
its rear end moves somewhat slower because of the lack 
of support by further strong energy release at the ori- 
gin. Therefore the fireball cone is preserved instead of 
being destroyed. However, preservation is a necessary 
condition for GRB production. The fireball can in fact 
be radially stretched before it reaches optically thin 
conditions. Assuming that the GRB duration, A^b, is 
defined by the time difference between the front and 
rear ends of the fireball reaching the "transparency ra- 
dius", e.g., 10^'^ cm, we can estimate, by extrapolating 
the results of our computed models, a burst duration 
Atb that might be significantly longer than the on-time 
(- 0.1 s) of the source: A^b = ^^-Stl^^ s. Note that this 
stretching can be the dominant contribution to Aib in 
case of source activity times of significantly less than 
one second as expected for the accretion timescale of 
post-merger disks. The latter timescale is expected to 
vary with the torus mass. It is set by the viscous trans- 
port within the compact accretion torus, in contrast to 
collapsars where the accretion disk is fed by an external 
reservoir of several solar masses of stellar matter and 
the accretion timescale is therefore determined by the 
collapse timescale of the massive, rotating star. Even 
if the accretion phase of the remnant BHs of NS+NS 
or NS+BH mergers lasts only fractions of a second, 
our simulations suggest that such events can account 
for the measured durations of short GRBs. Of course, 
our estimation of Aib can only be considered as an ex- 
ercise for demonstrating a fundamental possibility. It 
has no predictive power for GRB durations, because 
we need to extrapolate our hydrodynamic results over 
several orders of magnitude in radius. This is rather 
uncertain due to difficulties in extracting accurate val- 
ues for the fluid velocities from the computations in 
case of very high Lorentz factors and due to the ge- 
ometrical dilution of the spatial grid resolution with 
increasing radial distance. Extrapolation of our results 
from about 0.5 s 0.4 s) after the onset (shutdown) of 
the energy release to more than 300 s later also ignores 
how the fireball properties continue to change and how 
the long-time propagation and expansion of the fireball 
may depend on the structure of the ambient medium of 
the merger site. Moreover, the GRB emission might be 
shorter than our estimate if it is produced in a region 
that is smaller than the whole fireball. Nonthermal 



emission of radiation requires the dominant energy of 
the flow to be kinetic energy of relativistic baryons but 
not internal energy. Unfortunately our simulations had 
to be stopped before definite statements about the ter- 
minal fireball structure and the final ratio of its kinetic 
to internal energy were possible. 

— The total energy (internal plus kinetic) of the rel- 
ativistic fireball scales roughly linearly with the de- 
posited energy and saturates after the energy release 
has ended. A burst-like, short initial phase of en- 
ergy deposition, followed by a long-time, gradual de- 
cay turned out to channel somewhat less energy into 
the relativistic outflow and to be slightly less efficient 
in converting internal energy to kinetic energy than a 
constant rate of energy release with a sudden end. In 
both cases a minor part (a few per cent at most) of the 
energy was in form of kinetic energy after 0.5 seconds 
of computed evolution, with the tendency to continue 
rising. 

— The relativistic flow in type-B models reveals a highly 
inhomogeneous structure in both the radial and lat- 
eral directions. A small negative radial gradient of the 
radiation-dominated pressure corresponds to a radial 
increase of the Lorentz factor also in the shocked re- 
gion of the flow, where the average behavior of differ- 
ent quantities can roughly be described by power laws 
of the radius (though different from those of the ana- 
lytic model of .Levinson fc Eichler 2000). While pres- 
sure and Lorentz factor reveal only moderate varia- 
tions (factors of a few) , the density shows huge fluctu- 
ations up to three orders of magnitude. These inhomo- 
geneities are a relic of Kelvin-Helmholtz instabilities 
which are triggered by the interaction of the fireball 
with the torus medium. In lateral direction the struc- 
ture cannot be fitted by simple top- hat, power-law or 
Gaussian profiles but can posses off-axis maxima of the 
Lorentz factor. This feature results from a genuinely 
relativistic effect that occurs in the presence of tan- 
gential velocity components at a propagating shock. 
It is favored by the fact that the surfaces of constant 
energy deposition deviate from spheres leading to rel- 
atively high energy deposition rates at large radii off 
the symmetry axis. A relativistic shock with large tan- 
gential velocities causes the formation of a rarefaction 
wave that moves slowly from the lateral edges towards 
the axis of the polar outflow and allows the Lorentz 
factor in some of our B-modcls to become about 20% 
larger in the rarefied region than near the symmetry 
axis of the radial flow. At the same time the ultrarela- 
tivistic core, where Lorentz factors of several hundred 
are reached, is kept coUimated by a subsonic negative 
lateral velocity component while only the surrounding, 
mildly relativistic ejecta exhibit a very slow sideways 
expansion. Thus we expect the structure of the fire- 
ball to be essentially preserved during the subsequent 
evolution. 

— The lateral profile of the radially averaged energy of 
type-B models shows a central hot core (of w 3° — 9° 
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width) of nearly constant value, surrounded by a layer 
where the radially averaged energy drops proportional 
to . This profile arises naturally in all our mod- 
els, and agrees well with that of the so-called struc- 
tured jet model that has been used to explain the 
phenomenological correlati on between the GR B peak 
luminosity a nd pulse lags lIRorris et al.ll2000lj) of long 
GRBs (e.g., ISalmonsonll2n0f]l) or t he nearly constant 
energy inferred from observations bv^ Frail ^ 1^ ^iJ ( 2^^) 
and fro m theoretical models by lR,ossi et al.l (20021) 
and Zh ang fc MeszarosI l)2002l) . The lateral profile of 
the Lorentz factor displays a more complex structure 
whose functional form cannot be properly fitted by a 
9^^ or Gaussian decay (which is required in the stan- 
dard structured jet model) mainly because of off-axis 
maxima of the Lorentz factor. In case these profiles are 
able to survive until the fireball reaches transparency, 
the angular energetic profiles of short bursts may be 
very similar to those of the long bursts, but the Lorentz 
factor distribution might be more complex. 

Despite the coarsening of the computational grid with 
increasing radius, the fireball could be well followed from 
the shutdown of the energy source at tga (or during the 
long-time, gradual decay of the energy release, respec- 
tively) to the end of our simulations (at 0.5 s or 5isa) 
when the fireball reaches a radial extension of about 
1.5 X 10^° cm. For discussing gross observational features 
we therefore feel safe to assume that global fireball proper- 
ties can be sufficiently reliably extrapolated to much larger 
distances. Inhomogeneities on angular scales 9ih > F"^ 
should persist and are expected to be smoothed out only 
after th e Lorentz factor has dropped below (e.g., 
iKumar 'fc Granot 2003). The maximum Lorentz factors in 
the front of the outflow are higher than 500. In some of the 
type-B models regions in the fireball with F > 1000 are 
found at the end of our simulations. The terminal value of 
the Lorentz factor should depend on the ambient density 
of the merger site. Our results are obtained for a fireball 
that accelerates into a r"^ "* density profile (chosen mainly 
to overcome numerical problems in the construction of the 
initial configuration). However, the real situation may be 
more complex, e.g., a steeper density gradient might con- 
nect the torus "surface" (or a small cloud of baryons that 
surrounds the BH-torus system as the relic of the pre- 
ceding, violent merger history) with the much lower and 
almost constant density of the interstellar medium. 

With a half-opening angle of Oj ~ 5°-10° (15°-25°) 
the coUimated ultrarelativistic outflows of type-B mod- 
els having F > 100 (F > 10) at t = 0.5 seconds after the 
onset of the energy release by the GRB engine, cover 
a fraction = 1 - cos6'j ~ 0.4%-1.5% (< 10%) of 
the sky^. Assuming equal detectability from all directions 
within the opening angle this fraction implies about 100 

^ The numbers in brackets correspond to outflow with 
Lorentz factors 10 at 0.5 s, which shows still ongoing accel- 
eration so that much larger Lorentz factors can finally result, 
cf. Table H 



(more than 10)^ times more events than measured gamma- 
ray bursts. A rate of short GRBs of about 100 per year 
therefore requires an event rate per galaxy and year of 
10-^(/n/0.01)-i(iVg/109)-i where iVg is the number of 
visible galaxies. Comparing with estimated NS-I-NS and 
NS-I-BH merger rates, which are typically around 10~^ 
per year and gal axy (with about a factor of 10 or mor e 
uncertainty; e.g., iKaloeera et"al]l2004t IfYver et alJll999|) . 
we therefore conclude that a significant fraction of such 
mergers but probably not all, should produce GRB viable 
outflows. 

The edges of the ultrarelativistic (F > 100) jet core 
are very sharp in terms of the isotropic equivalent energy 
Eiso- Maximum (terminal) Lorentz factors of the order 
of 1000 suggest the potential to account for hard GRB 
spectra. Our simulations actually showed inhomogeneous 
and anisotropic, coUimated outflows with lateral variation 
of the Lorentz factor and of the apparent isotropic en- 
ergy. Hence, we do not expect equal observability from all 
positions in the beam direction. The maximum values of 
the apparent isotropic energy, E'iso, are found to be up 
to about 10*^^ erg at angles ^ 10° around the symmetry 
axis of the ultrarelativistic outflow, declining towards the 
outer wings of less relativistic ejecta. These numbers are 
obtained for an energy deposition rate of a few 10^° erg 
over a period of typically 100 ms. These are reasonable 
and not extreme values in view of model calculations for 
the energy release by neutrino- antineut rino- annihilation 
in case of post-rnergirig; BH accretion l|Ruffert fc Jankal 
ll999 H.Ta,nka et ablll 99flH S etiawan et a,].1 l2004h. Provided a, 
major fraction of the energy of the ultrarelativistic fireball 
gets converted to gamma-rays, our maximum isotropic 
equivalent energies are in good agreement with estimates 
based on a comparison of the energetics of short and long 
GRBs, sugge sting an approxima te fiuence-duration pro- 
portionality i Balazs et al.ll200^ . Since long bursts last 
typic ally about 50-100 times longer, a similar luminosity 
re.g.. lMao et 10119941: Lf^°'^ ~ L^°^^ - lO^i-s^ergs"!) 
implies an apparent energy which is around 10^^ erg for 
short bursts instead of ^ 10^'^ erg for long ones ijFrail et alJ 
Il997l) . 

Due to the lateral structure of the outfiow we also ex- 
pect some degree of variability between observed bursts, 
depending on the viewing angle relative to the system axis. 
Such differences will be superimposed on variations due 
to intrinsic properties of the binary systems and remnant 
BH-torus systems, e.g., associated with different masses 
and spins of the merging neutron stars or black holes. 
This finding should be taken into account in studies of 
the div ersity of short gamma-ray burst s like the recent 
one by ^Ross wog fc R.amirez-R.uia ()200,3l) . These authors 
also employed the assumption that the ultrarelativistic 
outflo w is confined as suggested by iLevinson fc Eichleil 
(l2000h . Our models, however, show a much different hy- 
drodynamic scenario for the fireball evolution and coUi- 
mation, in which the outfiow-torus interaction, relativis- 
tic shock effects, and Kelvin-Helmholtz instabilities play 
a crucial role. We therefore emphasize that conclusions 
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Fig. 31. Time evolution of the ratio of kinetic to total (internal plus kinetic) energy (left panel) of relativistic outflow 
in model BOl with Lorentz factor above 2, 10, 50 and 100, from top to bottom, respectively. The right panel shows 
the corresponding evolution of the total energy (solid lines) and of the kinetic energy (crosses) in the same order from 
top to bottom. 



based on grounds of the simplified picture developed by 
iLevinson fc Eichle^ l|2000() should be drawn with caution. 
Finally, we repeat that it is impossible to estimate the 
energetics of ultrarelativistic outflows without performing 
simulations that follow the complex hydrodynamics phe- 
nomena which develop in response to the deposition of en- 
ergy in the vicinity of the BH-torus system. A static anal- 
ysis of time slices for mass dist ribution and energy depo- 
sition by z/P-annihilation (e.g., iRosswog fc R.amirez-RAii3 
can therefore be misleading, in particular with re- 
spect to the outflow energetics and asymptotic Lorentz 
factor which are crucial for judging the viability of vD- 
annihilation for powering GRBs. 

Our simulations are only a preliminary step towards 
fully self-consistent models. So far we have investigated 
the relativistic hydrodynamic flow that is triggered by the 
deposition of energy near the BH-torus system. However, 
this energy release was prescribed according to a defined 
functional behavior instead of being linked to the neu- 
trino emission of the evolving accretion torus. We are 
in the process of removing this limitation by preparing 
two-dimensional relativistic simulations of the viscosity- 
driven evolution of BH-torus systems with a simplified, 
but consistent treatment of neutrino transport and vD- 
annihilation. 
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